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ABSTRACT : 

A new radio direction-finding (DF) technique which applies 
matched-filter theory to the DF problem is presented. The new 
technique, called Matched-Filter Doppler Direction Finding (MFD DF) 
is based on the doppler DF principle in which a rotating antenna, 
or its equivalent, superimposes periodic frequency modulation on the 
received carrier. The phase of the induced modulation contains 
information on the angle of arrival. It is shown that a bandpass 
filtering operation, which satisfies the phase-matching requirement 
of a matched filter, converts the frequency-modulated signal to an 
amplitude-modulated signal whose envelope is a periodic narrow pulse 
with no sidelobes. The relative time of occurrence of the pulse 
is a measure of angle of arrival. 

Included is the description of and results obtained with an 
experimental system used to confirm the analytical results, to study 
the effects of noise, and to explore alternatives available in the 
design of an operating system. Encouraging results were obtained 
using this same experimental system with a conventional operating 
doppler DF set. 

This task was supported by: Naval Ship Systems Command (Code 6050) 
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ABSTRACT 



A new radio direction-finding (DF) technique which 
applies matched-filter theory to the DP problem Is pre- 
sented. The new technique, called Matched-Filter Doppler 
Direction Finding (MFD DP) is based on the doppler DF 
principle in which a rotating antenna, or Its equivalent, 
superimposes periodic frequency modulation on the received 
carrier. The phase of the Induced modulation contains 
information on the angle of arrival. It Is shown that a 
bandpass filtering operation, which satisfies the phase- 
matching requirement of a matched filter, converts the 
frequency-modulated signal to an amplitude-modulated signal 
whose envelope is a periodic narrow pulse with no side- 
lobes. The relative time of occurrence of the pulse Is a 
measure of angle of arrival. 

Included is the description of and results obtained with 
an experimental system used to confirm the analytical 
results, to study the effects of noise, and to explore alter- 
natives available In the design of an operating system. 
Encouraging results were obtained using this same experi- 
mental system with a conventional operating doppler DF set. 
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I. INTRODUCTION 



This report is concerned with radio direction finding 
(DF), namely the ability to determine the direction of 
arrival of a radio wave. The direction of arrival includes 
azimuth angle (0-360 degrees) and elevation angle ( 0-90 
degrees); however, this paper deals almost exclusively 
with the azimuth angle of arrival (AOA) . In addition, an 
underlying assumption is that the receiver is far enough 
from the transmitter so that the arriving v/avefront is 
essentially planar. 

A. BACKGROUND 

As early as 1899 it v/as recognized that electromagnetic 
waves could be aimed or beamed in a specific direction by 
the use of antennas such as loops or spaced dipoles. For 
receiving purposes this meant that the direction of arrival 
of a radio wave could be determined by moving or steering 
a directive antenna beam in azimuth until a maximum signal 
was obtained. The beam could be steered by physically 
moving an antenna or by appropriately summing the outputs 
of the individual elements of an array. 

The science of radio direction finding has been the 
subject of much study since the early 1900 's as reflected, 
for example, in the extensive bibliography by Travers and 
Hlxon [1] which contains over 5000 entries. Among the 
Inventors of some of the first radio direction finders were 
such pioneers as Marconi and de Forrest [2]. 
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A major Improvement over early DF systems came wltl the 
development in the 1920-19^0 era of the Adcock DF system 
employing four or more spaced elements with an overall 
aperture of one-half wavelength or less. The Adcock system 
achieves an Instrumental accuracy of about five degrees [ 31 . 

Since the end of World War II there has been a trend 
toward larger aperture systems. For example, dopplei* DF 
systems, in which frequency modulation (FM) by a sinusoid 
is induced on the received signal by a rotating antenna or 
its equivalent, usually range in size from one-half to 
five wavelengths and are capable of accuracies of two 
degrees or less [3]* 

The Wullenweber DF system which was developed in Germany 
in the early 19^0's consists of a large circular array on 
the order of 10 wavelengths in diameter [4] . At any Instant 
of time only one segment of the array is used, and with 
proper phasing of the elements in use, this segment is made 
equivalent to a planar array. By commutation the planar 
array is in essence caused to move in a circle. An accuracy 
of one-half degree or less is possible with this system [3]* 

At microwave frequencies, antennas having aperture sizes 
of several wavelengths and with extremely directive patterns 
can be physically rotated to give accurate AOA information. 
An alternate method which is very accurate for large values 
of signal-to-nolse ratio (SNR) is the amplitude-sensing 
monopulse technique in which two identical directive 
antennas are placed some distance apart [ 51 . The pattern 
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produced by subtracting the output of one antenna from the 
other has a sharp null when the pair is aimed at the arriv- 
ing signal. 

B. CONTENTS AND SUMMARY 

In this report the DF problem is considered as a signal 
processing problem. In any scanning DF system an ideal 
waveform for determining AOA is a narrow pulse whose loca- 
tion along some time (azimuth) scale varies directly with 
the AOA. Figure 1 is an example of the processed output 
of the DF system considered in this report. In that figure 
the horizontal axis is calibrated to show AOA. It should 
be noted that the accuracy with which the location of the 




FIG. 1. PHOTOGRAPH OF THE TYPICAL OUTPUT 
OF A MFD DF SYSTEM 
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peak of any pulse can be detei*mlned is limited only by 
noise. In the presence of noise, Skolnik [6, p. 476-478] 
shows that measurement accuracy is directly proportional 
to the SNR and inversely proportional to the width of the 
pulse. From the well known Inverse relationship between 
a time function and its frequency representation as 
expressed in the Fourier transform, it is clear that the 
narrower a pulse in the time domain the larger its Land- 
width [7]. This basic principle, v/hlch underlies the 
theory of pulse compression, is applied here to the DF 
prob lem. 

To generate a narrow (short duration) pulse by signal 
processing, a useful spreading of the spectrum of the 
received signal is necessary. One way to broaden the 
signal spectrum is to pass the signal through a nonlinea^^ 
device. This is the basis for the DF technlqu'^ suggested 
by Terukhimovlch [8] in which the summed output of two 
antennas is applied to an nth-order nonlinear device to 
create higher hai-monlcs. The harmonics are amplified 
individually and then summed to produce a desired periodic 
train of narroxv radio-frequency pulses. This partlci lar 
technique has apparently received little attention. 

Another way of spreading the spectrum of an incominm 
signal is through use of the doppler principle. If a radio 
wave is received by an antenna element which rotates at a 
l'ix''d rate about a central point, and if the antenna xC 
omnidirectional in the plane of rotation, tlie Instantaneous 
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frequency of the received signal varies in a sinusoidal 
manner. It can be recognized from the doppler principle 
that the maximum positive frequency deviation occurs v/hen 
the antenna is moving perpendicularly toward the incoming 
wavefront; thus, there is a direct relationship bet^^reen 
the position of the antenna and the time of maximum instan- 
taneous frequency. The periodic PM signal produced in this 
manner not only has a spectrum whose bandwidth is propor- 
tional to the radius of rotation, but the phase of each of 
its discrete spectral components is directly related to AOA. 

Since a large bandwidth is only a necessary but not a 
sufficient condition for the generation of a short-duration 
pulse, other processing may be required. Processing the 
output of an antenna to improve its directivity or some 
other parameter is not a new concept. In fact, Ksienskl 
[9] surveys a number of signal-processing antennas. Hov;- 
ever, a DP system is specialized in that there is one 
particular desired output — an accurate indication of the 
AOA. Por this reason the design of a DP system may have 
fewer restrictions than the design of a directional 
receiving system. That is, a DP system need not Improve 
gain nor optimize any other parameter of reception; it must 
only provide accurate AOA Information. Por examiple, a null 
1 n an antenna pattern can be used to obtain an accurate 
bearing, but it is of no use in signal reception. 
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A principal result of this research Is the development 
of a Matched-Filter Doppler DF (MFD DF) technique which Is 
a simple means of processing the output of a rotating 
antenna to obtain the AOA of a radio wave. An MFD DF 
system, shown In general form In Fig. 2, extracts the AOA 
Information from the periodic FM signal by a linear filter- 
ing operation which yields a periodic amplitude-modulated 
signal consisting of a single recurring pulse with no 
sldelobes as shown In Fig. 1. The position of the voltage 
pulse relative to a reference time scale gives the AOA of 
the signal. Because the portion of the spectrum which is 
filtered satisfies the phase-matching requirement of a 
matched filter the new technique is called the Matched- 
Filter Doppler DF technique. 

Consideration is given to the use of frequency multi- 
plication which In the absence of noise Improves the 
ability to measure AOA accurately. The MFD DF concept Is 
simulated on a digital computer for modulation Indexes 
which correspond to a range of aperture diameters from 
eight to 1000 wa^'-elengths . 

To demonstrate the feasibility of the MFD DF principle 
and to investigate design alternatives in developing a 
practical system, an analog simulation was conducted using 
conventional electronic hardware. The same hardware is 
used in conjunction with a conventional doppler DF system 
to confirm the basic characteristics of the MFD DF tech- 
nique under operating conditions. 
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FIG. 2. BLOCK DlAGiRM OF A MFD DF SYSTEM 



c. 



PLAN OF THE REPORT 



Chapter II reviews the familiar doppler DF principle 
[ 10] [II] [12] , and develops in detail the theory of the new 
MFD DF teclinique. 

In order to apply time-domain signal-processing Tech- 
niques, such as matched filtering, to obtain spatial 
information, such as direction of arrival, the receiving 
antenna must convert spatial Information to temporal form. 

As shown in Chapter III this is exactly what occurs when 
a linear antenna (or array) is rotated about a central 
point. In addition, this chapter unifies a number of 
previously known results of antenna theory, applies them 
to rotating antennas, and, for the first time, extends the 
results to rotating signal-processing DF antennas in 
general, and the MFD DF system in particular. 

Chapter IV demonstrates the physical realizability of 
the MFD DF scheme by presenting results obtained wit), a 
hardware simulation system. Physical limitations are 
considered as well as the effect of noise. The last section 
in this chapter gives the results of an experiment in which 
the simulation hardware is interfaced with a conventional 
doppler DF set. 

Chapter V is a theoretical study of the effect on a 
MFD DF system produced by the simultaneous reception of two 
signals of the same frequency at different azimuth angles. 



Chapter VI is a summary of the report and lists 



recommendations for future work. Appendices A, B, 
C support results presented in the report. A list 
references is provided. 



and 

of 
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II. THE MATCHED-FILTER DOPPLER DF TECHNIQUE 



This chapter reviews the doppler DF principle and then 
introduces the matched-filter doppler DF technique in which 
the AOA of a radio wave is extracted from a frequency- 
modulated carrier by a linear filtering process. The tech- 
nique relies on the fact that m.odulation is imposed on the 
received radio wave by the motion of the receiving element. 
The periodic frequency-modulated signal at the antenna out- 
put is converted to a periodic amplitude-modulated signal 
consisting of a narrow voltage pulse with no sldelobeo. 

The position of the pulse relative to a fixed reference 
time is directly related to the AOA. 

A. THE DOPPLER DF PRINCIPLE 

Consider a radio wave of carrier frequency f^ received 
by an antenna element which rotates at a constant rate 
f^ about a central point and v/hlch is omnidirectional in 
the plane of rotation as shown in Fig. 3- For simplicity 
the radio wave is assumed to lie in the plane of rotation. 
If the azimuthal AOA is 6^^ the output of the antenna is 
given by [Ref. 13] 

s(t) = Ccos[oj t + B(slnoj t - 9 )] (1) 

c o 

where C is a constant, to = 2iTf = 2tt/A , to = 27rf and the 

c c c r r 

modulation index or peak phase deviation is 

B = 2tt£/A^ (2) 
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Reference 




FIG. 3. PLAN VIEW OF A ROTATING 

antenna element 



where I is the radius of the circular antenna path.^ The 
same doppler effect can be achieved with a circular array 
of receiving elements samipled sequentially. 

Equation (1) Is In the familiar form of frequency or 
phase modulation by a sinusoid. As the antenna rotates 
the Instantaneous frequency of the antenna output varies 
due to the doppler effect. As the receiving element travels 
towards and then away from the Incoming wavefront the 
received signal Is frequency modulated In a sinusoidal 
fashion. This doppler effect can be viewed as phase 

^If the signal is not traveling in the horizontal plane 
of rotation and arrives at an angle of elevation a the only 
effect of interest here is to make the modulation Index 
3 = 27rA( cosa)/A^ . 
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modulation by noting that as the antenna travels around 
the circle, the path length between the transmitter and 
the receiving element changes. This causes the phase of 
the received signal to vary sinusoidally relative to the 
phase of the signal at the center of the circular path. 

The peak phase deviation 6 Is therefore quite naturally 
equal to the radius of rotation In radians as given by (2). 
More Importantly, for DF purposes, regardless of vdiether 
the antenna output Is considered frequency modulation or 
phase modulation, the relative phase of the modulating 
sinusoid Is directly related to the AOA of the signal. 

This doppler effect Is the basis for conventional doppler 
DF systems, and It Is also the basis for the matched- 
filter doppler DF technique presented in this report. 

B. CONVEIITIONAL DOPPLER DF SYSTEMS 

Conventional doppler DF systems were first proposed 
In the late 19^0's by Budenbom [10] and Boulet [11] in 
this country and by Earp and Godfrey [12] In Great Britain. 

The block diagram of a typical conventional doppler DF 
system Is given In Fig. The output of a rotating antenna 

or an equivalent commutated array Is frequency demodulated 
yielding a sinusoidal voltage 

s^(t) “ 3co^cos(oj^t - e^) . (3) 

The phase of s^(t) is compared with the phase of a like- 
frequency reference sinusoid which Is associated v.'lth a 
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FIG. 4. BIOCK DIAGRAim OF A CONVhNTIONAL DOPPLER DF SYSIH 



fixed direction, typically true north. As shown in Fig. 4, 
the dc voltage level out of the phase comparator is a linear 
measure of the AOA. Conventional doppler DF sets are 
usually of medium-aperture size (circular path diameter of 
one-half to five wavelengths), and their performance In 
terms of accuracy Is reported to be two degrees or better 
[3]. Conventional doppler systems cannot measure the AOA 
of weak signals in the presence of strong signals or noise 
due to the capture effect or small-signal suppression effect 
of FM receivers [^] [1^] . 

The matched-filter doppler DF technique presented in 
the next section requires no frequency or phase demodulators 
(nonlinear devices). The MFD DF system is completely 
linear, easily realizable, and from Figs. 1 and 2 It is 
seen that the MFD DF system provides AOA information In a 
simple and useful form suitable for direct display. 

C. THE MATCHED-FILTER DOPPLER DF SYSTEM 

In a scanning DF system, a repetitive narrow pulse 
whose occurrence time (phase) varies linearly with AOA is 
a useful waveform. In this section matched-filter process- 
ing Is considered as a means of obtaining such a pulse 
from a carrier frequency modulated by a sinusoid. 

1 . Matched-Filter Concepts 

It is useful here to review the properties of 
matched filters [151 • If a function v(t) has a Fourier 
transform V(f) given by 
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oo 



( 4 ) 



V(f) = / v(t)e = |V(f) 

— OO 

where 4'(f) is the argument of the complex frequency function, 
then the transfer function of a filter matched i-o v(t) is 
given by 



K^|V(f) 






g-j2TTfX 



(5) 



where is a constant and t is a fixed time delay exper- 
ienced by the signal v/hen passed through the filter. The 
matched filter is also called a conjugate filter because, 
except for the delay and amplitude factors, the filter 
transfer function is the conjugate of V(f); therefore, a 
matched filter: 

(1) Scales the magnitude of each spectral component 
of v(t) by the factor |V(f)| (amplitude matching) 
and, 

(2) Places all of the spectral components of v(t) in 
phase at t = T due to the conjugate property (phase 
matching) . 

To use these fundamental properties of matched 
filters to provide DF information, the spectral representa- 
tion of FM by a sinusoid Is now considered. 

2 . Spectral Representation of FM by a Sinusoid 

The function s(t) given in (1) Is periodic with 
period 2iT/m^ and hence can be expanded In a Fourier series 
as [16, p. 125-130] , 
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s(t) = J (3)cos(jot + Z J, (3)cos[(o) +kaj)t-k6] 

O C ^ C O 

oo , 

+ Z (-1)^ J, (3)cos [(U)^ -ku )t + ke^] (6) 
k=l ^ c r o 

where is the Bessel function of the first kind of 

order k and argument 3. In (6) and equations following, C 
has been set equal to one. It is assumed throughout the 
following development that f^<<f^. 

Equation (6) shows that FM by a sinusoid can be repre- 
sented by a carrier-frequency term and an infinite numiber 
of side-frequency terms, each separated from its neighbor 
by a frequency increment equal to f^ Hz. Corresponding 
upper and lower side frequencies have the same magnitude. 
Bessel functions also have the property that beyond the 
value of order k equal to the argument 3, the magnitude of 
the Bessel function and hence of the corresponding frequency 
component decreases rapidly and asymptotically to zero [17] . 
This means that although the series is infinite, there are 
a finite number of significant frequency components. 

When considered as phasors, the frequency components 
given in (6) have relative phase-angle relationships which 
vary with time; however, there are specific instants. of 
time when all of the phasors have either a zero or l 80 -degree 
relative phase angle. These Instants are Important here 
because they are directly related to the AOA of the received 
signal. For example. Fig. 5 shows the magnitudes of the 
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FIG. 5. SPECTRUM OF FM BY A SINUSOID WHEN 
Reference time t =0 /r . 



coefficients of the Fourier series components of s(t) when 
6 = 25 . (A value of 6 of 25 corresponds In Fig. 3 to a 
radius of about four wavelengths.) The sign of each 
coefficient is indicated by its projection above (positive) 
or below (negative) the frequency axis. Conventionally, 
such a diagram as Fig. 5 is called the spectrum of the 
signal . 

Because the series which represents s(t) consists 
solely of cosine functions and since all of the coefficients 
(Bessel functions) are real. Fig. 5 can be used to cor»vey 
information about the relative phase angles between the 
various frequency components (phasors). However, since 
each phasor rotates at a ’lifferent rate. Fig. 5 only gives 
valid phase-angle relatlor ships when the principal value of 
the argument of e ve ry cosine function in (6) is the same. 

The first instant at which this happens is t = 0 if the 
cosine terms have no initial phase angle (0^ = 0). Since 
s(t) has a phase lag of 0. radians, the instant at which 
the relative phase angle r lationships of Fig. 5 exist is 



t = 0 /u . At t 
o o r ( 



o 



0 /m 
o 



s(0 /to ) = J (3)costo (0 /to ) + 
o r o c o r 



k=l 



E J, ( 3) costo ( 0 /to ) 
, k ^ cor 



00 



(-1)^ Jj^(3)costo^(0^/(o^) 



+ E 
k=l 



(7) 



and all of the spectral components (phasors) with positive 
coefficients (above the frequency axis) have a zero-degree 



phac'=' angle relative to each other, while those spectral 
cciT.ponents with negative coefficients (below the frequency 
axis) have a l80-degree relative phase angle. As seen 
from the second term of (6) and (7), the spectral com.ponents 
of frequency greater than the carrier frequency have a 
positive or negative coefficient according to the sign of 
t"he Bessel function with which each is associated. When 
3 = 25, all Bessel functions of order 20 or more are posi- 
tive; hence, all frequency components of value f^ + 20f^ 
or more are in relative phase conjunction at t =6 . 

Because s(t) has a period of 2 tt/o)^, the same relative 
phase-angle relationships shown in Fig. 5 recur at instants 



t 



m 



6 

o 



+ 2mTT 
0 ) 



r 



m=0,l, ...n. (8) 



The preceding analysis was concerned primarily with 
the characteristics of the frequency components above the 
carrier frequency. The reason for this is now given. Each 
side frequency below the carrier is equal in magnitude to 
its counterpart above the carrier and differs only in 
relative phase angle at t = because of the (-1)'^ 

factor in the third term of (6) and (7). However, if a new 
reference time t' is selected where 



t ' = t - tt/o) , 

which corresponds to a displacement of one-half of a rota- 
tion cycle of the antenna of Fig. 8, then the antenna output 
can be expressed as 
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s ( t ) = J (3)cosw ( t ) + E (-1) J, (6)cos[co ( t ) +kco t-kfi ] 

0 ) o c 0 ) , k c oo r o 

r r k=l r 



+ E J, (3)cos[o) (t ^)-koa t + k6 ] 

k=l ^ ‘ 0 )^ r o 



(9) 



The spectrum of (9) is given in Fig. 6 and it reflects 

the change in location of the (-1)^ factor from the third 

term in (6) to the second term in (9)- By similarity with 

the argument presented for s(t) as represented in Fig. 5j 

Fig. 6 gives the relative phase angles (zero or l80 degrees) 

between the spectral components of s(t') at time t^ = 0^/o)^ . 

At t', which corresponds to t = t + tt/o) 
o^ ^ o r 

0+TT 0+7T 00 , 0+7T 

) = J (3)coso) (— ) + E (-1)^ J, (3)cosoo (— ^ ) 

oj o c oo ,, k C O) 

r r k=l r 

00 0 + IT 

+ E J, (3)cosoo^(^- — ) (10) 

k=l r 

It is clear that the relative phase angle relat^ionshij.'s 
among the phasors (frequency components) above the carrier 
frequency at t^(Fig. 5)> are transposed to the frequency 
components below the carrier frequency at t^ + tt/oo^ (Fig. 6). 
The significance of this result is that there are only two 
Instants during each rotation or scanning cycle and 

(0^+tt)/oo^ during the first cycle) when all of the spectral 
components of the antenna output have either a zero or a 
l80-degree relative phase angle; these two instants are 
directly related to the AOA of the signal. 
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FTG. 6. SPECTRUM OF FM BY A SINUSOID WHEN 



3 . Determining Angle of Arrival 



Having reviewed matched-filter theory and investi- 
gated the spectral characteristics of FM by a sinusoid, 
this section describes a realizable filter which is phase 
matched to part of the spectrum of the antenna output 
voltage. The filter output can be used directly -- no 
further processing required — to obtain the AOA of a radio 
wave . 

Consider first the output of a filter matched to 
s(t) of (1). The matched-filter output is easy to calculate 
mathematically by multiplying each coefficient of the 
Fourier series expansion of s(t) given in (6) by Itself. 

This operation yields a function which satisfies both the 
amplitude-matching and phase-matching requirements of a 
matched filter. The computed envelope of one cycle of the 
filter output is shown in Fig. 7-a for B =25. This envelope 
can be used to determine AOA directly since it consists of 
a narrow main lobe whose position along the time axis 
depends on the AOA of the signal. That is, the time axis 
can be graduated in degrees of azimuth for DF purposes. 

The voltage of Fig. 7-a repeats once for eacli rotation cycle 
of the antenna. It is interesting to note that the 3-db 
width of the main lobe is 5.04 degrees (considering one 
scanning cycle equivalent to 360 degrees), which is nearly 
two degrees narrower than the calculated main lobe of an 
antenna pattern produced by a uniformly excited linear 
antenna of the same aperture. 
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(a) Amplitude and phase matching 




(b) Phase matching only 

FIG. 7. ENVELOPE OF ONE CYCLE OF THE PERIODIC OUTPUT 
OF A FILTER MATCHED TO A CARRIER FM BY A SINUSOID . 

3 is 25. 
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It is observed from Fig. 5 that various adjacent 
frequency components of s(t) have phase angles which differ 
by IT radians; therefore any filter matched to s(t) must 
also Introduce ir-radian phase changes in a similar manner 
due to the phase-matching property. The realization of such 
a filter is impractical. Phase matching is in genera] more 
important in the generation of narrow pulses than is ampli- 
tude matching. This is shown in Fig. 7-b which displays 
the calculated envelope of one cycle of the output of a 
phase-matched rectangular bandpass filter for B = 25. For 
DF purposes the waveform of Fig. 7-b is just as useful as 
that of Fig. 7-a. However, in both cases the presentation 
may be undesirable because of the number and magnitude of 
sidelobes present. 

Difficulty in realizing a filter matched to the 
spectrum of Fig. 5 and the sldelobe structure of the matched- 
filter output (Fig. 7 ) prompts consideration of compromise 
solutions. When B > 10 the spectrum of FM by a sinusoid 
always contains one large contiguous group of components 
(phasors) which have the same relative phase angle at the 
Instant t^ = When B = 25 these like-phase components 

appear as the group shown isolated in Fig. 8. A portion of 
the spectrum, therefore, already satisfies the phase-matching 
requirement of a matched-filter output. If these components 
are removed as a group from the rest of the spectrum, the 
resulting waveform should display some of the properties 
frequently seen in matched-filter waveforms -- a large 
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IN-PHASr GROUP OF INTEREST 



amplitude at one instant of time with relatively smaller 
amplitudes at all other times during each cycle of antenna 
rotation. The In-phase group having the largest number of 
significant components, which for 3 = 25 consist of the 
20th through 31st harmonics of f^ above f^, can be selected 
with a bandpass filter whose output is given from (6) as 

31 2 

s,(t) = E J, ( 3 )cos [ (w +kw )t - k6 ] . (11) 

1 ^ c r o 

By the use of trigonometric identities s^(t) can 
be expressed as 

s^(t) = E(t)cos [o)^t + (J)(t) ] (12) 

where E(t) is the amplitude-modulation factor (envelope) 
and (J)(t) is a phase modulation term (see Appendix A for 
details ) . 

The waveform E(t) of Interest, shown in Fig. ?-a for 
one scanning cycle, not only has a main peak corresponding 
in time to the AOA of the signal (l80° in this case), but 
it is also a smooth pulse free of sidelobes. As the trans- 
mitter azimuth changes the pulse position moves along the 
time axis. Figure 9-b is a photograph of the undetected 
output of an experimental system used to simulate the MFD DF 
system. Experimental results are presented in detail in 
Chapter IV. 



2 

In this and other examples the last significant compo- 
nent is the last component that has any discernable effect 
on the calculated waveform s^^(t). 



H2 




l.o: G.DJ 



(a) Computed envelope. 




(b) Undetected output of experimental system. 

Two scanning cycles are displayed. 

FIG. 9. ENVELOPE OF THE ISOLATED IN-PHASE GROUP OF 
FREQUENCY CCMPONENTS SHOWN IN FIG. 8 
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Because of the spectral symmetry discussed In 



section 2 and shown In Fig. 6, there Is a corresponding 
group of frequencies below f^ which are also In phase 
relative to each other once per cycle but at a later time 
corresponding to an additional half-cycle of antenna 
rotation. Obviously this lower-frequency group can be 
bandpass filtered to produce a waveform Identical to that 
shown In Fig. 9, except for the fixed time delay of 
seconds. Computations and experimental results verify 
this conclusion. 

It Is Instructive to consider the In-phase spectral 
components as a group of sinusoids which are added. 
Gulllemln [l8] calls waveforms such as s^(t) In (11) fre- 
quency groups. These particular groups have uniformly 
spaced discrete-frequency components and all frequency 
components are In phase at some reference time. Gulllemln 
explores some of the properties of frequency groups, par- 
ticularly the envelopes of groups with equal-amplitude 
components because they can be expressed mathematically In 
closed form with easily predictable behavior. For example, 
a group of 12 sinusoids of equal amplitude separated by 
Af H and with a median frequency f much greater than Af 
can be expressed as 



P(t) 



. rl2(27TAf).. , 

sin [ -t ] 

! r2lTAf^ 1 
Sin [ — t ] 



cos2irf t 
c 



(13) 



2 



The envelope of p(t) as shown In Fig. 10 has a prin- 
cipal maximum of width l/(6Af) v;hlch repeats every l/(Af) 
seconds. Minima occur every l/(12Af) seconds except when 
this corresponds to a multiple of l/(Af). Secondary "axlma 
(siaelobes) appear about halfway between the secondary minima. 
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FIG. 10. COMPUTED ENVELOPE OF AN EQUAL -AMPLITUDE, 

IN-PHASE FREQUENCY GROUP 

Unfortunately, groups such as those encountered in 
the spectrum of Fig. 5 which have members of unequal magni- 
tudes can not be so simply expressed as (13), and their 
envelopes are not as readily predicted from their series 
form (11) or closed form (12). For example, as sh^'iwn 
previous l.y, the group of 12 Iri-phase spectral components 
associated with FM by a sinusoid when 3 is 25 (Fig. 8) is 
expressed as the series of (11). 



^5 



The envelope E(t) of this group as given in (12) 



has a maximum at multiples of 1/f^, but other specific 
details are not obvious unless the envelope is plotted. 
Despite the differences between equal-amplitude frequency 
groups and those of Interest here, Guillemln's results can 
be used to obtain Insight into general envelope properties 
of frequency groups. 



insight into the behavior of a frequency group. In phasor 
form, a group such as s^(t) of (11) adds together in the 
complex plane to form a resultant phasor of magnitude E(t) 
and phase angle 4>(t). To illustrate this. Pig. 11 shows 
a typical group of five phasors at two Instants of time. 



The phasor diagram is also a useful tool to gain 




FIG. 11. PHASOR DIAGRAMS 



Because s^(t) is an In-phase group and because, the 
frequency components which form s^(t) differ from each 
other by f^ Hz, once and only once every 1/fj, seconds all 
of the phasors add constructively to form the large pulse 
sliown in Fig. 9 . It is interesting to note that there are 
no times when smaller groups among the frequency components 
of s^(t) constructively Interfere to form sidelobes. 

In summary, it has been shown that a rotating 
antenna element superimposes sinusoidal FM on a received 
carrier. The simple process of bandpass filtering a group 
of in-phase spectral components of the FM signal produces a 
periodic AM waveform with no dlscernable sidelobes. The 
peak of the lobe occurs at a time linearly related to the 
AOA of the signal. Since the in-phase group of frequencies 
satisfies the phase-matching property of a matched-filter 
output, the new technique is called Matched-Filter Doppler 
Direction Finding (MFD DF). These theoretical results are 
supported by results obtained with an experimental MFD DF 
system. These experimental results are presented in 
Chapter IV. 

k . Dependence of DF Accuracy on System Parameters 

In the absence of noise the time of occurrence of 
the peak of a pulse can be measured in theory with unlimited 
accuracy. However, the ease with which this value can be 
ascertained, particularly by an operator, is dependent on 
the shape of the pulse. Furthermore, Skolnik [6, p. 
shows that the measurement error in the presence of noise 



(at a high SNR) is inversely proportional to the square 
root of the SNR and directly proportional to the half-power 
width of the pulse. In the limit as the SNR becomes very 
large, the error becomes negligible. This section explores 
the properties of FM by a sinusoid and the matched-filter 
technique, particularly the effect of those parameters 
which affect the shape of the output pulse of a MFD DF 
system . 

The two parameters of s(t) of (1) that can be con- 
trolled at the receiving site are scanning rate f^ and 
modulation index 6. As will be shown later, the value of 
f^ affects system design, but its only effect on the 
spectrum of s(t) is to change the separation between the 
frequency components . The value of modulation index g 
given by (2) has a marked effect on the spectrum. As 6 
Increases the total number of significant sidebands 
Increases. For 6 greater than 50, approximately 2gsig- 
nlficant sidebands exist [171. Figure 12 shows the spectrum 
of s(t) for values of g of 50 and 100 (only components 
above f^ are shown) . The characteristic of primary impor- 
tance here is that the number of components in the last 
group of in-phase members and hence the bandwidth of that 
group increases with g. It is also interesting that the 
difference in magnitudes of adjacent components in that 
group decreases with g. 

As expected from Guillemln's results and as shown 
in Fig. 13, obtained by evaluating equations similar to 
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FIG. 12. SPECTRA OF FM BY A SINUSOID 
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(a) 6 = 25, Half- power width = 43.5° 




(b) 6 = 50, Half- power width = 33.8° 



FIG, 13. COMPUTED ENVELOPES OF THE OUTPUT OF AN IDEAL 
MFD DF SYSTEM AS B INCREASES 



50 




(c) 6 = 100, Half- power width = 27.5° 




(d) B = 200, Half- power width = 21.8° 
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(e) 6 = 400, Half-power width = 17.20^ 




(f) S e 800, Half-power width = 13,8° 
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(g) 6 = 1600, Half- power width = 11,5° 




(h) B = 3200, Half- power width = 8.2° 
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(11) with the aid of a digital computer, the main lobe of 
the envelopes of the in-phase frequency groups narrows with 
increasing 3 and no dlscernable sldelobes appear. Experi- 
mental results confirm this as seen in Pig. 14 which shows 




FIG. 14. OUTPUT OF THE EXPERIMENTAL MFD DF SYSTEM 
WHEN 0 = 100 

the output of an experimental MFD DF system when 3 is 100 
(Chapter IV gives details). Figure 15 summarizes these 
results by simultaneously displaying the 3-db widths of 
the main lobe of the envelopes as well as the number of 
members in the in-phase group as a function of 3 • It is 
concluded then that as 3 increases, the width of the output 
pulse of the MFD DF system narrows which is deslreable for 
measurement purposes because the time at which the pulse 
occurs is more dlscernable to an operator. 
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5 . Increased B by Frequency Multiplication 



For a given carrier wavelength B is linearly 

related to the radius of the path of the rotating antenna 
or equivalent array [Eq. (2)]. Practical limits are soon 
reached in attempting to increase the value of B by 
increasing the array diameter. However, frequency multi- 
plication [l6, p. 132 - 13 ^] provides a means of increasing 
B without increasing the size of the antenna aperture. 

For a carrier frequency modulated by a sinusoid given by 

s(t) = Ccos[co^t + Bsin ( w^t-e^ ) ] , (1) 

frequency multiplication by a factor n produces 

Sn(t) = Ceos [nm^t + nB(slno)^t - 6^)1 • (1^) 

The carrier frequency is increased by the factor n, but the 
spectral components of s^(t) are still separated by f^ Hz. 
More importantly, the modulation index is now nB and the 
relative phase of the modulating signal, which corresponds 
to the AOA, is preserved. Frequency multiplication for 
the application of interest has the same effect as increas- 
ing the size of the receiving aperture. 

The effect of noise on the technique of frequency 
multiplication will be covered in Chapter IV; however, even 
in the absense of noise. Fig. 15 shows that the envelope 
3-dB width decreases at a rather slow rate with increasing 
B so that practically there is undoubtedly a limited number 
of useful stages of frequency multiplication. 
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FIG. 15. GRAPH OF THE 3-dB PULSE WIDTH AND THE NUMBER 
OF SIGNIFICANT COMPONENTS IN THE IN-PHASE GROUP VS 6 



There is no theoretical limit to the value of B 
possible and hence to the width of the output pulse (after 
bandpass filtering the in-phase group). This can be inter- 
preted, for DF purposes at least, as corresponding to 
superdirectivity in the time domain. The term super- 
directivity is used in the sense that, given a signal from 
a finite-aperture antenna, a waveform can in theory be 
produced which has a main pulse of duration less than the 
duration of a pulse formed by the scanning beamwldth of 
the antenna pattern produced by a uniformly excited aper- 
ture of the same size. This corresponds to the broad 
definition of superdlrectlvlty given by Tucker [191 . The 
concept of superdirectivity will be considered in Chapter 
III with regard to the processing of the output of rotating 
antennas . 
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III. RADIO DIRECTION FINDING WITH ROTATING ANTENNAS 



A single rotating receiving element, used to this point 
to explain the MFD DF technique. Is but one of many possible 
forms of rotating antennas. This chapter considers rotating 
antennas In general, with emphasis on their direction- 
finding capabilities. The purpose of the chapter Is 

1. To unify a number of previously known results of 
antenna theory and to apply these results to rotating 
antennas . 

2. To extend these results to Include signal processing 
of the output of rotating antennas, and 

3. To apply these results to the specific case of 
MFD DF. 

Antennas which rotate at a constant rate convert spatial 
Information to temporal Information. Once this conversion 
has occurred, signal processing techniques can be applied 
to produce a waveform from which desired spatial Informa- 
tion, such as the bearing of a radio source, can be accur- 
ately extracted. In the analysis of rotating antennas It 
proves useful to consider a Fourier transform representation 
of the azimuthal distribution of far- field sources. In the 
transform domain, antennas are shown to act as "filters" 
to the far- field sources. 

The general theory of rotating antennas Is applied to 
the case of a single rotating element and It Is argued that 
the output of such an antenna, by virtue of Its spectral 
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characteristics, contains as much directional information 
as the pattern of a highly directive array of the same 
aperture size. This is confirmed by the results of 
Chapter II. 

A. THEORY OP ROTATING ANTENNAS 
1 . Basic Formulation 

Consider a one-dimensional aperture or an antenna 
array with a far-fleld pattern y(9). For a receiving 
antenna this means that its output v(t) due to a point 
source radiating a sinusoid of frequency f^ and located at 
an azimuth 6^ is given in phasor notation by 

TT 

v(t) = / y(0) (1>(0 - e^)e d0 (15-a) 

-TT 

where 6(0 -9^) is ths ordinary Dirac delta function. 

Prom the property of delta functions, (15-a) becomes 

j (“pt) 

v(t) = y(0^)e (15-b) 

Although frequently not considered, y(0) may have a phase 
characteristic as well as an amplitude characteristic. We 
observe from (15-b) that v experiences a sinusoidal time 
variation. In the following development, this time varia- 
tion is implied when writing v as a function of azimuthal 

r 

angle . 

If the antenna is pivoted at its center and may be 
rotated, an additional variable, angle y of* rotation, is 



59 



Is Introduced and the antenna pattern as a function of two 
variables becomes y(0-Y). 

In the general case, which Includes distributed 
far-field sources given by z(0), (15-a) becomes 

v(y) = / y(e-Y)z(0)d0 (16) 

-7T 

where v(y) and z(y) are understood to be sinusoidal time 
functions [ 20 ] . 

If y(0) is assumed to be an even function, (16) can 
be written as 

v(y) = / y(Y-0)z(0)d0 (17) 

-TT 

The right-hand member is in the form of a convolu- 
tion Integral except for the finite limits of integration. 

In (17) v(y) is a periodic function because y(Y-6) is 
periodic in y and v(y) can therefore be expanded in a 
Fourier series which has a Fourier transform consisting of 
discrete components. However, for the following development 
it is useful to define a modified function, v^(y) having a 
continuous Fourier transform. The results are then related 
to the periodic case. To establish v^(y) the function 
y^(0) is Introduced such that 

f y ( 0 ) -■n'<.0£ + ir 

y]^(0) = < o (18) 

0 elsewhere . 

A single exception is the case where y(0) is constant. 
In this case y^(0) = y(6). 
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Equation (17) can then be rewritten with extended 



limits of Integration as 



oo 



v^(y) = f y^(Y - 0) z( 9)d0 . 



(19) 



This formulation Is similar though somewhat more general 
than that given by Bracewell and Roberts [20] . Since the 
modified variables y^(6) and v^(y) as well as z(6) have 
continuous Fourier transforms. It proves useful to examine 
their Fourier transform representation. 

2 . Fourier Transform Representation 



By analogy to the Fourier transform of a time function, ? 
may be considered a frequency variable. It has no physical 
significance comparable with the physical meaning associated 
with frequencies arising from the Fourier transform of a 
time function. 

The variable ^ is related to spatial frequency n 
which arises when the Fourier transform of an antenna 
pattern, expressed as a function of u = sln6, is taken as 



Consider the Fourier transform of a function of 6 



(or y) such as z(9) 



Z(C) = / z(6)e . 



( 20 ) 



G(n) = g(u)e . 



( 21 ) 



_00 



6l 



It must be emphasized that the variables C and n are 
not linearly related except when 6=sln0. Therefore, In the 
antenna literature, results given in terms of spatial fre- 
quency n are not, in general, directly applicable to the 
^-frequency representation [21] . One well-known result con- 
cerning spatial frequency comes from the fact that the far- 
field pattern g(u) of a one-dlmenslonal unidirectional antenna 
oriented along the x-axls Is proportional to the Fourier 
transform of the aperture excitation function D(x/A), 

g(u) a r d(x/A) (22) 

— 00 

where 6 Is measured from a perpendicular to the x-axls as 
shown in Fig. l6 [22]. 




-L/2 T 



FIG. 16. ONE -DIMENS ZONAL APERTURE AND RELATED GEOMETRY 
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The Inverse transform is given by 



D(x/X) « / du (23) 

— oo 

Because this integral includes a region of u which 
is physically meaningless, that is |u|>l, the following 
result is meaningful only for directive antennas for which 
g(u)-^0 when | y | ^1 . 

By noting the similarity between (21) and (23) it 
is seen that directive antennas have spatial frequency 
functions G(ri) which are proportional to the spatial 
excitation D(x/X) of the aperture. Because the aperture 
is absolutely limited in extent, the spatial frequencies 
are absolutely limited also [21] . It must be emphasized 
that this result applies to the ^-frequency representation 
only in cases where g(u) approaches zero while sin 9 and 9 
are Interchangeable. 

Having defined the Fourier transform of a function 
of azimuth angle, (19) can be considered again. The 
Fourier transform of v^(9) in (19) is given from the con- 
volution theorem by 

V^(C) = • Z(?) (24) 

If the antenna excitation is considered a linear process 
then the antenna itself may be considered a "filter" Y^(?) 
acting upon the ^-spectral characteristics of the far field, 
Z(^). The "filter" characteristics can be controlled to 
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some extent by the designer, so it is Important to determine 
what characteristics are desirable in a DF system. 

In order to obtain an exact azimuth bearing and 
perfect resolution, an ideal antenna pattern is a delta 
function 



y^(e-Y) = 6(6 -y) • (25) 

In the C domain this gives 

Y^(^) = (26-a) 

from which 

|y2(c)| = 1 • (26-b) 

Correspondingly, a point source at azimuth 6^ is 
described by 

z(e) = 6(0-9^) (27) 

and in the x, domain, 

-J27t6 

Z(^) = e ° (28-a) 

and 

1Z(^)| = 1 . (28-b) 

The ideal point source therefore has c-f'i’squency 
components extending over the entire x, spectrum and the 
transfer function of the ideal DF antenna preserves all 
these components. This results in a delta function for 
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thG output v^(6) which gives the azimuth angle of the sounce 
precisely . 

Consider by means of an example what occurs In the 
case of a real antenna. Suppose the antenna pattern Is 
given by 



y^Ce) 



g-e^/2a^ 

/Stto 



(29) 



where a Is such that y^(6) Is 
The transform of (29) 



very small when 
Is given by 



I 6 I 2.TT/2 . 



Y^(?) = Ke 



2 2 2 
-2t\ a 



(30) 



If we consider for simplicity a point source at zero degrees 
of azimuth then. 



z(6) = 6(6) 



(31) 



The antenna output is given In the C domain by 

V^(C) = Y^(S) • Z(i;) , (32-a) 

2 2 2 

V^(0 = Ke"^"" ^ • 1 . (32-b) 

Comparing this result with the Ideal case It can 
be seen that a real antenna has a "filter" characteristic 
which attenuates C-spectral components of the far-fleld 
source In proportion to C* Obviously, as the antenna 
pattern y^(0) gets narrower (a decreases) the corresponding 
transfer function Y^(c) becomes broader. Conventionally, 
the directivity of an antenna pattern Is Increased by making 
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the antenna larger in aperture which has, of course, 
practical limitations. Since real antennas have patterns 
which extend over a finite range of 6 their transfer 
functions in the S domain are never all-pass as in the 
ideal case but taper off in magnitude with ? in some 
fashion, depending on the specific antenna. This contrasts 
with the spatial frequency approach in which the spatial 
frequency representation of an antenna pattern is absolutly 
limited along the q-axis. The example just given in 
(29)-(32) shows that the ^-frequency transform of an antenna 
pattern does not have absolute limits though its magnitude 
does decrease directly with c;. Davies and Longstaff [231 
arrived at this same basic result for the case of a rotating 
linear array. 

As mentioned previously, an antenna pattern ideal 
for DF is one with an arbitrarily narrow beamwldth. The 
quest for unlimited directivity leads to consideration of 
superdirective antennas which are broadly defined as anten- 
nas which exceed the directivity of a uniformly excited 
antenna of the same size [19]* The theory of superdirective 
(or supergain) antennas shows that there is no limit to 
the directivity that can be achieved with a finite aperture 
antenna; however, supergain is difficult to achieve in 
practice [2^1 ] . 

Kslenskl [21] has shown that superdirective antennas 
cannot be characterized in terms of spatial frequency u 
because the simple Fourier transform relationship between 
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the antenna excitation function and the far field pattern 
given by (22) and (23) loses Its physical connotation. That 
ISj classical superdlrectlve antennas have far-fleld 
patterns which have finite theoretical values, outside of 
the range of real angles (|u|>l). Hence, the assumption 
which allowed use of limits of ±°° In ( 23 ) does not hold. 
On the other hand Davies and Longstaff [231 showed by means 
of an example which is directly applicable to Z frequencies 
that conventional techniques for designing supergain arrays 
have the effect of attenuating the lower ? components to 
make them comparable in amplitude to the higher ? components 
or In effect broadening the "filter" passband by emphasizing 
the higher C frequencies. Consequently, this suggests an 
alternate definition of supergain as any technique which 
Increases the directivity of an antenna beyond that of an 
equivalent uniformly excited antenna by Increasing Its 
C-spectrum passband without increasing its aperture. 

3 . Time-Domain Processing 

In the preceding paragraphs antennas have been con- 
sidered In terms of their C-frequency characteristics. If 
it were possible to convert K frequencies to real frequen- 
cies, the converted output of an antenna could be processed 
to achieve an overall transfer function including the 
antenna and the processor which produces in the time domain 
a waveform with good directional characteristics , perhaps 
even superdlrectlvlty . Assuming that the conversion to the 
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real frequency domain Is possible, as will be demonstrated, 
such a system is shown In block diagram form as Pig. 17- 




Antenna Converter Processor 



FIG. 17. BLOCK DIAGRAM OF A FREQUENCY-PROCESSING 
RECEIVING SYSTEM 

From basic transform algebra, for a linear system, 

W(f) = H(f) Y^(f) Z(f) (33) 

As discussed previously. It is desirable that W(f) 
equal Z(f), so that w(t) is a delta function. This requires 
that 

1 = H(f) Y^(f) (34-a) 

from which 

H(f) = 1/Y^(f) . (3^-b) 

This makes H(f) the Inverse filter of Y^(f). In a similar 
context, Davies and Longstaff [23] refer to the Inverse 
filter implicitly In proposing a filter that equalizes all 
frequency components to the same magnitude. For the example 
given previously by (29)-(32), 
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Therefore, from (3^-b) 



H(f) 



2tt a f 



K 



( 36 ) 



Ii(f) is a perfectly acceptable processor transfer function 
in a noiseless environment and the processor output w(t) cor- 
responds to an arbitrarily narrow voltage pulse; however, 
in the presence of white noise, the inverse filter for 
functions such as Y^(f) accentuates the noise spectrum 
indefinitely at high frequencies and is thus not practical. 
Qualitatively, it should be noted that the Inverse filter 
broadens the input spectrum by accentuating the higher 
frequencies and hence satisfies the alternate definition 
of supergain proposed in section A. 2 of this chapter. 

The best known compromise to the inverse filter is 
the matched filter which for a given signal maximizes the 
SNR at the filter output at some Instant of time based on 
knowledge of the noise spectral density [15]. Usually, 
additive white gausslan noise is assumed. 

It is necessary at this point to digress briefly 
to relate the results obtained using y^^CS), a truncated 
version of the actual antenna pattern y(6), to the actual 
situation. By using y^(6) it was possible to make the 
preceding ^-frequency development by means of continuous 
Fourier transforms and basic transfer-function algebra. 



since y(y) and the antenna output v(y) are periodic 
functions of the rotation angle y> they can be expanded In 
a Fourier series having a spectrum of discrete components 
In the C-frequency domain. The envelope of these discrete 
spectral components corresponds to the hypothetical 
V^(?) or ¥2^(5), just as the envelope of the Fourier trans- 
form of a rectangular pulse train corresponds to the trans- 
form of a single pulse. With this distinction In mond 
the conclusions arrived at In the preceding sections are 
applicable . 

^ . Conversion of g frequencies to Real Frequencies 

The process of converting c frequencies to real 
frequencies Is straightforward. When a linear antenna or 
array rotates at a constant rate its output is given 

from (17) by 

v(w^t) = / y(to^t-e)z(e)de (37) 

-TT 

What previously was a function of y becomes a 
function of time which repeats every 2 ir/u)^ seconds. The 
Fourier transform of the output Is no longer In the 
^-frequency domain but rather In the real-frequency domain. 
Uniform rotation of the antenna can be thought of as pro- 
ducing a transformation of variable from C to f. As 
discussed previously the antenna output can then be pro- 
cessed as desired. 



The repetition period Is only 7r/o3 seconds for linear 
ariays which arc usually bldlrectlonal^ln their pattern. 
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B. ONE ROTATING ELEMENT 



This section applies the general theory of rotating 
antennas to the specific case of a single rotating element 
which is omnidirectional in the plane of propagation. It 
is shown that the output of such an antenna contains 
directional information which, of course, confirms the 
results obtained in Chapter II using a different approach. 

A single, fixed vertical dipole has a horizontal antenna 
pattern given by 

y^ce) = K (38) 

In the 5 domain this becomes 

y^(0 = K 6(0 (39) 

As a filter to the ^-spectrum of a far-field point source 

this antenna has an extremely narrow bandpass. Since 

directive antennas have been associated with broad 

^-spectral bandpasses, one would expect the fixed dipole 

on this basis to have little directivity in the horizontal 

plane. The output of the dipole due to a point source at 

6 = 0 is, in phasor notation 

o 

j U) t 

v(t) = Ke ^ (40) 

As expected, neither the envelope nor the phase of the 
output voltage contains information about the signal AOA. 

However, if a dipole which is omnidirectional in the 
plane of the circle of Fig. l8 is placed a distance I from 
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FIG. 18. PLAN VIEW OF THE ROTATABLE ELEMENT 



a central reference point and oriented initially in azimuth 
as shown, its excitation pattern is given in phasor notation 
as 



1 

y(e) 




2-irZ 

X 



c 



sin 9 



(^ 1 ) 



where the phase is measured relative to that of the voltage 
measured at the center of rotation in Fig. 18. If the 
element is displaced y radians counter-clockwise around 
the circle, the pattern becom.es 

J-j^sin(y-e) 

y(e-y) = Ke ^ . (il2) 



If the element rotates at a constant rate f^, the 
pattern becomes a function of time and azimuth. 
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J^^sln(w^t-e) 



y(e-w^t) = Ke 



c 



(^ 3 ) 



The Fourier transform of y(0,to^t) is 



y(f) = K I 

k=-cx> 



)e J^®6(f-kf^) . 



( 33 ) 



The bandwidth of Y(f) Is proportional to the size of the 



Since, from Fourier transform theory, pulse duration Is pro- 
portional to the Inverse of bandwidth, then the minimum- 
pulse duration obtainable from the response of this antenna 



to a far-fleld point source Is approximately seconds. 

The pulse repeats every 1/fj, seconds, which corresponds to 
one 360 -degree scan In azimuth, so the minimum width of 
the pulse In radians Is given by 



As pointed out by Davies and Longstaff, who were considering 
the outer elements of a rotating array, this corresponds 
roughly to the common rule of thumb for the pattern of a 
directive linear array, namely, that the beamwidth In 
radians is approximately equal to the wavelength divided 
by the array length, 21, 




TriTf” A 
r c 




beamwidth 



wavelength _ _c 



array length 2H 
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This suggests that the output of a single rotating 
element contains enough directional information to deter- 
mine the AOA with at least as much accuracy as a directional 
array of comparable aperture. 

This fact is confirmed theoretically in Fig. 7-a which 
shows that the output waveform of a filter matched to the 
output of a rotating antenna with a diameter of about eight 
wavelengths (3 = 25) has a 3 dB width of 5.0 degrees. The 
calculated beamwldth of a uniformly excited linear array 
eight wavelengths long is 7.2 degrees. 
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IV. HARDWARE SIMULATION AND PERFORMANCE OF A MATCHED- 

FILTER DOPPLER DF SYSTEM 



This chapter reports on the results of a hardware 
simulation of the MFD DF system. The experimental system 
also defines some of the limitations of the new technique 
arising from, the use of conventional electronic hardware. 

A series of experiments produced data on the effect of 
noise on the system for both multiplicative and non- 
multlplicatlve processing. The experimental system was 
also used with an operational doppler DF system. The 
results of this test confirm the principal features of 
the matched-filter technique. 

A. SYSTEM DESCRIPTION 

The experimental system is shown in block diagram form 
in Fig. 19 . The output of a rotating antenna (or equiva- 
lent array) is simulated by a voltage controlled oscillator 
(VCO) which is frequency modulated at a rate of f^ Hz 
corresponding to the scanning rate of an actual doppler DF 
system. The processing begins with an optional frequency- 
multiplier section capable of multiplication by a factor 
of two or four (Fig. 20). A balanced mixer and a variable 
local oscillator are used to shift the signal spectrum in 
order to bring any desired group of frequencies within the 
filter passband. 

The filter output, which contains the desired AM, is 
available for display Immediately, or if desired the RF 
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FIG. 19. BLOCK DIAGR;«^ OF THE EXPERIMENTAL 
MFD DF SYSTEM 




a 

c 
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can be first passed to a simple diode envelope detector 
to obtain the AM voltage. 

The critical hardware component is the bandpass filter. 
The filter must pass only the desired spectral components, 
and the filter phase characteristic must be relatively 
linear. The required filter bandwidth is determined by 
the separation of f^ Hz of the discrete frequency compo- 
nents and by the number of significant components in the 
in-phase group. This number of components depends on the 
modulation index (B) of the signal applied to the filter. 

A modulation or scanning rate of f^ = 100 Hz is used because 
it is comparable to the scanning rates of some conventional 
doppler DF systems. 

After limited experimentation with tuned, ceramic, and 
mechanical filters, a Clevite ceramic filter (TL-10D9-20A) 
with a nominal passband of 10 kHz centered at H55 kHz was 
found to give results which most nearly approach the wave- 
forms obtained by digital computation. The characteristics 
of the filter used are shown in Fig. 21. For some experi- 
ments, a Clevite filter (TL-4D8A) with a narrower passband 
of about ^ kHz was also used. The characteristics of the 
4-kHz filter are shown in Fig. 22. 

B. CONFIRMATION OF THEORY 

1 . Processing Without Frequency Multiplication 

The important characteristics of a MFD DF system 
are demonstrated by the following example. The VCO of 
Fig. 19 generates the following signal. 
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G. 21. AMPLITUDE AND PHASE -VS -FRFQUENC 
CHARACTERISTICS OF THE 10-kHz CERAMIC 
BANDPASS FILTER 



Attenuation ( dB ) 




Frequency, kHz 

FIG. 22. J^MPLITUDE AND PHASE -VS -FREQUENCY 
CHARACTERISTICS OF THE 4 -kHz CERAMIC 
BANDPASS FILTER 
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Relative Phase Lag , Degrees 



(^ 5 ) 



s(t) = Ccos((o t + gslnco t) 
c r 

where = 429 kHz and f = 100 Hz. 

c r 

By adjusting the carrier frequency the signal 
spectrum is shifted until the in-phase group of spectral 
components (Pig. 8) is in the filter passband (frequency 
multiplier bypassed) . The resulting filter output before 
and after envelope detection is shown in Pig. 23- The 
experimental and calculated envelopes are remarkably 
similar as can be seen by comparing Pig. 23 with the com- 
puted output shown in Pig. 24. The 3-Hb width of the 
envelope pulse is about 50 degrees in the experimental case 




FIG. 24. COMPUTED ENVELOPE OF THE RESPONSE OF AN 
IDEAL MFD DF SYSTEM 
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(a) Predetection envelope of bandpass filtered 
FM by a sinusoid. g is 25. 

X Scale: 2 ms/cm 




(b) Detected envelope of bandpass filtered 
FM by a sinusoid. 3 is 25. 

FIG. 23. PHOTOGRAPHS OF THE RESPONSE OF THE EXPERIMENTAL MFD DF SYSTEM 
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as compared with ^3 degrees for the calculated envelope 
pulse. The difference Is due to two factors: 

(1) In order to exclude the undeslred components some 
of the desired frequencies are supressed due to 
the non-ideal "skirt" of the filter. 

(2) The nonlinear phase characteristic of the ceramic 
filter tends to broaden the main lobe of the 
envelope. The narrowing of the main lobe with 
Increasing 6 Is seen In Pig. 25 which shows the 
experimental system output for B's of 50 and 100. 
These pulses can be compared with their calculated 
versions shown In Fig. 13-c and 13-d of Chapter II. 

2 . Spectral Symmetry 

It was previously shown that the effect of filter- 
ing that group of harmonics below the carrier frequency f^ 
which corresponds to the In-phase group above f^ simply 
shifts the pulse position In time an amount corresponding 
to a l80-degree shift In AOA. This Is clearly seen In 
Fig. 26 which shows the envelope of the In-phase group 
below the carrier frequency for the signal given In (45). 

The main pulse Is shifted one-half of a modulation cycle 
or 2.5 cm to the left of the position of the pulse shown 
In Fig. 23 . The distortion In Fig. 26 Is due to the 
asymmetry of the ceramic-filter characteristics, l.e., the 
lower cutoff (filter skirt) Is used to separate the In-phase 
components above the carrier frequency , while the upper 
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a. a is 50. 




b. 6 is 100. 



FIG. 25. PHOTOGRAPHS OF THE RESPONSE OF THE 
EXPERIMENTAL MFD DF SYSTEM AS 6 INCREASES 
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FIG. 26. PHOTOGRAPH OF THE RESPONSE OF THE EXPERIMENTAL MFD DF SYSTEM TO 
THE IN-PHASE GROUP OF FREQUENCY COMPONENTS BELOW THE CARRIER FREQUENCY 

.cutoff Isolates the corresponding group below the carrier 
frequency . 

3 . Processing with Frequency Multiplication 

The technique of frequency multiplication to 
increase 3 and thereby Improve bearing determination by 
producing a narrower pulse is dem.onstrated by starting with 
a signal having a 3 of 25, and frequency multiplying by 
a factor of four to increase 3 to 100 prior to filtering. 

The filtered RF voltage shown in Fig. 27 has a 3-db width 
of 31 degrees compared with 50 degrees when 3 = 25 (Fig. 

23), and 27.5 degrees in the ideal case (Fig. 13-c). As 
expected, the result is almost identical to that obtained 
by generating a signal with a 3 of 100 and filtering direct- 
ly with no multiplication as shown in Fig. 25-b. 
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(a) Predetection envelope of bandpass filtered FM 
by a sinusoid after multiplication by four. 

3 increased from 25 to 100. 




(b) Same as (a) after envelope detection. 

FIG. 27. PHOTOGRAPH OF THE RESPONSE OF THE EXPERIMENTAL MFD DF SYSTEM 
WHICH INCLUDES FREQUENCY MULTIPLICATION BY A FACTOR OF FOUR 
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C. IMPERFECT FILTERING EFFECTS 



In earlier sections of this report it has been shov/n 
that the last group of in-phase spectral components above 
(or below) the carrier frequency contains AOA information 
in its periodic envelope. This section examines the effect 
of Imperfect filtering of the in-phase group from the 
remaining components of the spectrum of the rotating 
antenna output. 

The inclusion in the filter passband of harmonics of 
lower frequency than the in-phase group is shown to produce 
the underlrable effect of splitting the main lobe of the 
output envelope and reducing its amplitude. It will also 
be seen that the exclusion of harmonics in the in-phase 
group reduces the amplitude of the main lobe and broadens 
it, but this exclusion does not alter the position of the 
main lobe. Hence no error is Introduced in the AOA 
measurement . 

Imperfect filtering may occur for the following reasons 

(1) Realizable bandpass filters do not have infinite 
cutoff characteristics. 

(2) Local oscillator mlstunlng can shift the spectrum 
presented to the bandpass filter. 

Regardless of the modulation index 6 the spectral com- 
ponents of a signal which is frequency modulated by a 
sinusoid possess a pattern of alternating in-phase and out- 
of-phase groups of side frequencies as shown earlier in 
Figs. 5 and 6. For experimental convenience, the 
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Investigation of the effect of Imperfect filtering Is 
carried out for the case of B = 118; but the results are 
generally applicable for all values of B greater than about 
10 . 

The partial spectrum for B = 118 Is shown In Fig. 28. 

The last In-phase group begins with the 109th harmonic of 
the scanning frequency f^ above the carrier frequency f^ 
and extends to the last significant harmonic which Is the 
127th. The calculated envelope of the In-phase group Is 
shown In Fig. 29* As harmonics below the 109th are Included 
In the filter passband they have an adverse effect on the 
envelope as seen In Fig. 30 because the first group of 
undeslred harmonics (IO3-IO8) are relatively out-of-phase 




FIG. 29. COMPUTED ENVELOPE OF THE RESPONSE OF AN 
IDEAL MFD DF SYSTEM 
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FIG. 28. PARTIAI. SPECTRUM OF FM BY A SINUSOID 



n 
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n — 




(b) Envelope when harmonics 105 - 108 are included 

FIG 30 THEORETICAL- RESPONSE OF A MFD DF SYSTEM AS 
UKTOESIRED FREQUENCY COMPONENTS ARE PROCESSED 
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(c) Envelope when harmonics 103 - 108 are included 




(d) Envelope when harmonics 101 - 108 are included 



FIG, 30, THEORETICAL RESPONSE OF A MFD DF SYSTEM AS 
UNDESIRED FREQUENCY COMPONENTS ARE PROCESSED 
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at the instant that the in-phase group is in phase. This 
reduces the peak value of the envelope and, in addition, 
as more harmonics are included there is an interaction 
between the components in the in-phase and out-of-phase 
groups which produces a splitting of the main lobe of the 
envelope . 

Because realizable filters have finite cutoff character- 
istics, it is impossible to duplicate experimentally the 
condition from which the computer-drawn envelopes of Fig. 

30 are obtained. It is possible, however, to shift the 
sinusoidal PM spectrum in the experimental system described 
earlier in such a way as to Include in the passband a number 
of spectral components below the 109th harmonic. The 
results are presented in Figs. 31-33 as undesired frequency 
components are brought into the filter passband. They 
substantially confirm the computer results. Figure 31-a 
shows an output closely resembling the theoretical result 
given as Fig. 29. Figures 32-a and 33-a show the emergence 
of the second peak which did not grow to the same size as 
the first, probably due to the nonlinearity of the phase 
characteristic of the 10-kHz ceramic filter. Along with 
each photograph is shown the spectrum which appears at the 
output of the filter (Figs. 31-b - 33-b). The dotted 
spectral lines indicate the magnitudes prior to filtering. 

Of as much Interest as the inclusion of undeslred 
harmonics is the exclusion of some of the desired ones. 

This is explored for the ideal case in the computer-drawn 
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Relative amplitude 




(a) Envelope of the spectrum shown in Fig. 31(b), 




(b) Spectrum of FM by a sinusoid after bandpass filtering. 
3 is 118. 

FIG. 31. RESPONSE OF THE EXPERIMENTAL MFD DF SYSTEM 
AS UNDESIRED FREQUENCY COMPONENTS ARE PROCESSED 
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Relative amplitude 




(a) Envelope of the spectrum shown in Fig. 32(b). 




f requency 



(b) Spectrum of FM by a sinusoid after bandpass filtering, 
g is 118. 

FIG. 32. RESPONSE OF THE EXPERIMENTAL MFD DF SYSTEM 
AS UNDESIRED FREQUENCY COMPONENTS ARE PROCESSED 
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Relative amplitude 




(a) Envelope of the spectrum shown in Fig. 33(b). 




frequency 



(b) Spectrum of FM by a sinusoid after bandpass filtering. 
6 is 118. 



FIG. 33. RESPONSE OF THE EXPERIMENTAL MFD DF SYSTEM 
AS UNDESIRED FREQUENCY COMPONENTS ARE PROCESSED 
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graphs shown in Fig. 3^ in which in-phase components start- 
ing with the 109 th harmonic are removed two at a time up 
to the ll4th harmonic. The results are more predictable 
because the remaining members still form an in-phase group 
and from Guillemin's frequency-group theory it can be 
anticipated that as the number of members decreases the 
principal lobe will broaden proportionately [I 8 ] . This 
general prediction is borne out by Fig. 3^ which also 
suggests that the absence of sldelobes of the envelope 
is caused by the tapering off in magnitude of J^(B) with n. 

Attenuation of the lower- frequency components in the 
in-phase group reduces the peak level of the envelope and 
increases the pulse width, but the envelope remains with- 
out sidelobes. These analytical results are substantially 
confirmed experimentally as shown in Figs. 35-37 which show 
the effect of a non-ideal bandpass filter when the FM 
spectrum is shifted to exclude or greatly attenuate, some 
of the desired in-phase components (109-127). 

D. NOISE EFFECTS 

This section explores the effects of bandllmited 
gaussian noise on the performance of the experimental MFD DF 
system. First, the system hardware is described. Next, 
matched- filter processing without frequency multiplication 
is considered, theoretically, for a slgnal-plus-nolse 
input, and an Improvement in SNR is predicted. Experimental 
results for this case are presented. In theory, multiplic- 
ative processing is shown to degrade SNR. It is estimated 
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'(a) Response when harmonics 111 through 129 are included. 




(b) Response when harmonics 113 through 126 are included. 

FIG. 34. THEORETICAL RESPONSE OF A MFD DF SYSTEM 

AS SCME DESIRED FREQUENCY COMPONENTS ARE NOT PROCESSED 
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Relative amplitude 




(a) Envelope of the spectrum shown in Fig. 35(b). 




f requency 



(b) Spectrum of FM by a sinusoid after bandpass filtering. 
3 is 118. 



FIG. 35. RESPONSE OF THE EXPERIMENTAL MFD DF SYSTEM AS 
DESIRED FREQUENCY COMPONENTS ARE NOT PROCESSED 
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Relative amplitude 




(a) Envelope of the spectrum shown in Fig. 36(b). 




(b) Spectrum of FM by a sinusoid after bandpass filtering, 
g is 118. 



FIG, 36. RESPONSE OF THE EXPERIMENTAL MFD DF SYSTEM AS 
DESIRED FREQUENCY COMPONENTS ARE NOT PROCESSED 
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(a) Envelope of the spectrum shown on Fig. 37(b). 




(b) Spectrum of FM by a sinusoid after bandpass filtering. 
0 is 118. 

FIG. 37. RESPONSE OF THE EXPERIMENTAL MFD DF SYSTEM AS 
DESIRED FREQUENCY COMPONENTS ARE NOT PROCESSED 
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that a 6-dB loss in SNR occurs each time the frequency 
Is doubled. Experimental results confirm this loss in SNR. 
These data are presented and the performance of a system 
which uses frequency multiplication is compared with one 
that does not when the input is signal plus noise. 

1 . Equipment 

The experimental system is shown in Fig. 38 . The 
system is basically the same as that used to demonstrate 
the MFD DF concept except for the addition of a noise 
generator and noise bandpass filters. 

The output of a wideband gaussian noise generator 
(Elgenco model 603A) is bandpass filtered to simulate the 
bandllmlted noise in a radio receiver. The VCO output and the 
noise are combined in a summing network and then applied 
to the processing hardware. 

When the frequency-multiplier section is bypassed, 
the VCO output is centered at 363 kHz, and the local 
oscillator is set about ^55 kHz higher. A center frequency 
of 363 kHz ensures that the original signal will not inter- 
fere with the difference frequency at the mixer output. 

The attenuation-vs-frequency characteristics of the circuit 
used to limit the bandwidth of the output of the noise 
source (centered at 363 kHz) is shown in Fig. 39-a- 

When the frequency multiplier is added, a VCO 
frequency of 455 kHz poses no problem and permits the use 
of a ceramic filter (centered at 455 kHz) as an optional 
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FIG. 38. BLOCK DIAGRAM OF THE EXPERIMENTAL MFD DF 
SYSTEM USED TO STUDY NOISE EFFECTS 




Attenuation (dB) Attenuation (dB) 
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Frequency, kHz 

FIG, 39. AMPLITUDE -VS -FREQUENCY CHARACTERISTICS OF 
BANDPASS NOISE FILTERS 
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means of obtaining bandlimlted noise. The tuned-circuit 
noise-filter characteristics (455 kHz center frequency) are 
given in Fig. 39-b . All signal and noise voltages were 
obtained using a Hewlett-Packard HP 3400A true-RMS voltmeter. 

2 . Processing Without Frequency Multiplication 
a. Theory 

Assuming that white, bandlimlted, gaussian noise 
is added to the signal, the theoretical SNR at the output 
of the system relative to the input SNR can be determined 
directly from the processing scheme. Assume an ideal filter 
passes a fraction d of the total bandwidth, and therefore 
of the total noise power, containing only the desired 
frequencies which constitute a fraction d' of the total 
signal power. The change in SNR caused by such filtering 
is given by 

= 10 log ^ (46) 

The value of increases slowly with 6 as shown in Fig. 40. 
This can be understood qualitatively since the total band- 
width of a signal frequency modulated by a sinusoid is 
approximately equal to 2 g, while the number of components 
in the desired in-phase group Increases much more slowly 
with g. 

Since the information on signal AOA is contained 
in the time of occurrence of the peak of the MFD DF system 
output, then a useful SNR is the ratio of peak signal 
voltage to rms noise voltage. This new ratio can be 
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FIG. 40. GRAPH OF THK SNR IMPROVEMENT 



obtained by defining a term which relates the peak 
envelope voltage to the rms envelope voltage: 

I 2 = 20 log (il 7 ) 

^rms 

Fig. ^1, obtained by computer analysis, shows 
that I 2 Increases as a function of 6. 

The terms and I 2 are Improvement terms and, 
as such, are Independent of the input SNR. Input SNR 
depends on receiver bandwidth. We assume a receiver band- 
width which passes all of the significant spectral com- 
ponents of the PM signal appearing at the output of the 
receiving antenna. Therefore the assumed input bandwidth 
is a function of the modulation index 6 and the scanning 
rate f^. For any given B, reduction of f^ decreases the 
spacing between spectral components thus reducing the over- 
all bandwidth. However, f^ cannot be reduced Indefinitely 
because of the need to filter a selected group of side 
frequencies. If f^ is too small, a desirable bandpass 
filter characteristic cannot be realized. The optimum 
bandwidth is determined therefore by engineering consider- 
ations, primarily those imposed by the construction state 
of the art of narrow bandpass filters having acceptable 
(linear) phase characteristics. 

In summary, processing without multiplication 
produces a theoretical Improvement in SNR which is directly 
related to the modulation index B. In turn, B is directly 
related to the radius of rotation (aperture) of the 
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GRAPH OF THE SNR IMPROVEMENT 



receiving antenna. System SNR varies Inversely with 
Since system accuracy Is proportional to SNR It Is 
desirable to maximize 3 (Increase aperture size) and to 
minimize consistent with practical considerations such 
as the availability of suitable bandpass filters, 
b. Experimental Results 

The effect of bandllmlted gausslan noise on 
the experimental MFD DF system which uses no frequency 
multiplication Is the subject of three experiments. The 
first experiment uses a FM signal with a B of 25 to which 
bandllmlted noise Is added. The mixer output Is passed to 
a ^-kHz bandpass filter (Clevlte TL 4D8-A). This filter 
extracts the desired 12 side frequencies occupying 1.2 kHz 
of bandwidth. The Input SNR Is varied from zero to 10 dB. 
The results are presented as indicated in Table I. The 
first photograph shows two cycles of the envelope of the 
output of the bandpass filter; the other photograph indi- 
cates this envelope voltage after linear-diode detection 
and low-pass filtering with an RC filter having a time 
constant (4.7 ms) comparable with the scanning period 
(10.0 ms). This large value of time constant provides 
integration which improves the ability of an operator to 
determine the time of occurrence of the pulse. The family 
of photographs associated with Table I is shown In Fig. 42 
Because the filter bandpass (4 kHz) is not exactly equal 
to the bandwidth of the in-phase group of side frequencies 
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TABLE I 



EFFECT OF NOISE ON THE PERFORMANCE OF 
AN EXPERIMENTAL MFD DF SYSTEM 

Basic Parameters : 

(a) Processing without frequency multiplication 

(b) 3 = 25 

(c) ^-kHz bandpass filter (TL ilD8-A) 



Figure 


INPUT 




OUTPUT 




Signal 

(dB) 


Noise 

(dB) 


SNR 


Signal 

(dB) 


Noise 

(dB) 


SNR 


42-a,b 


-13.0 


_ oo 


cx> 


-15.0 


_ OO 


OO 


42-c,d 


-13.0 


-13.0 


0 


-15.0 


-12.0 


-3.0 


^2-e,f 


-13.0 


-16.0 


3.0 


-15.0 


-14.2 


0.8 


^2-g,h 


-13.0 


-19.0 


6.0 


-15.0 


-17.0 


2.0 


^2-i 


-13.0 


-23.0 


10.0 


-15.0 


-21.5 


6.5 
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(a) 




(b) 

FIG. 42. PHOTOGRAPHS OF THE RESPONSE OF THE EXPERIMENTAL 
MFD DF SYSTEM FOR VARIOUS VALUES OF INPUT SNR 
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(d) 
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(g) 




(h) 
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there Is actually a decrease in SNR associated with the 
filtering operation. (The term of (^6) is negative.) 

For the second experiment the modulation index 
6 of the signal is increased to 100. All other sys^ ( m 
parameters are unchanged. The ^-kHz filter still passes 
the l8 in-phase side frequencies. Once again the input 
SNR is varied. The results are presented in Table II and 
Pig. ^3* The effect of the waveform Improvement term I 2 
due to the higher 6 can be seen by comparing any two photo- 
graphs in Pigs. ^2 and 43 which have approximately equal 
values of output SNR. For example, compare Fig. 42-g,h 
with Fig. 43 - 6 , f. Since the 4-kHz bandpass filter is much 
larger than necessary, the improvement I^ caused by filter- 
ing is about the same for the case 6 = 25 as for the case 
3 = 100; therefore the primary difference in display is 
due to the larger value of I 2 which is associated with 
larger values of 3. 

The use of a detector circuit with a large time 
constant broadens the peak of the envelope-detected pulse. 
This broadening destroys a portion of the accuracy improve- 
ment which results from an increased value of 3. With 
large values of input SNR such long integration timies are 
not desirable. 

The third experiment in this series is similar 
to the second with the exception that a 10-kHz bandpass 
filter ( Clevite TL 10D9-20A) replaces the 4-kHz filter. The 
results given in Table III and Pig. 44 provide both visual 
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TABLE II 



EFFECT OF NOISE ON THE PERFORMANCE OF 
AN EXPERIMENTAL MFD DF SYSTEM 



Basic Parameters : 

(a) Processing without frequency multiplication 

(b) 6 = 100 

(c) i^-kHz bandpass filter (TL 4d8-A) 



Figure 




INPUT 




OUTPUT 




Signal 

(dB) 


Noise 

(dB) 


SNR 


Signal 

(dB) 


Noise 

(dB) 


SNR 


i^3-a,b 


-13.0 


_ oo 


00 


-15.5 


— 00 


00 


H3-c,d 


-19.0 


-13.0 


-6 


-18.5 


-18.5 


0 


(U 

1 


-13.0 


-13.0 


0 


-15.5 


-18.7 


3.2 


^3-g,h 


-13.0 


-19.0 


6 


-15.5 


-23.7 


8.2 


^3-1. J 


-13.0 


-23.0 


10 


-15.5 


-28.5 


13.0 
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(a) 




FIG. 43. PHOTOGRAPHS OF THE RESPONSE OF THE EXPERIMENTAL 
MFD DF SYSTEM FOR VARIOUS VALUES OF INPUT SNR 
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(c) 




(d) 
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(e) 




(f) 
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(h) 
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TABLE III 



EFFECT OF NOISE ON THE PERFORMANCE OF 
AN EXPERIMENTAL MFD DF SYSTEM 

Basic Parameters : 

(a) Processing without frequency multiplication 

(b) 6 = 100 

(c) 10-kHz bandpass filter (TL10D9-20A) 



Figure 


INPUT 




OUTPUT 




Signal 

(dB) 


Noise 

(dB) 


SNR 


Signal 

(dB) 


Noise 

(dB) 


SNR 


^4-a ,b 


- 13.0 


— 00 


00 


-15 


— 00 


00 


44-c ,d 


- 19.0 


- 13.0 


- 6.0 


- 21.0 


-lh.5 


- 5.5 


44-e ,f 


- 13.0 


- 13.0 


0 


- 15.0 


- 14.5 


- 0.5 


4^-g,h 


- 13.0 


- 19.0 


6 . 0 


- 15.0 


- 20.0 


5.0 
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(a) 




(b) 



FIG. 44. PHOTOGRAPHS OF THE RESPONSE OF THE EXPERIMENTAL 
MFD DF SYSTEM FOR VARIOUS VALUES OF INPUT SNR 
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(d) 
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(f) 
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and experimental evidence of the decrease in fNH and conse- 
quent degradation in system accuracy caused by the additional 
noise power which passes through the increased filter 
bandpass. The 10-kHz passband contains about 2h times more 
noise power than the 4-kHz passband. The effect of this 
additional noise is seen by comparing any two cases from 
Tables II and III (Figs. 43 and 44) having the same input 
SNR. In general there is a 3 to 4 dB decrease in the out- 
put SNR when the 10-kHz filter is used. 

3 . Processing With Frequency Multiplication 

Frequency multiplication artificially Increases 
the aperture size of a MFD DF system by Increasing the 
modulation index 3 of the FM signal produced by the 
antenna rotation. This section explores the effect of 
additive, bandllmlted gausslan noise when frequency 
multiplication is employed, 
a. Theory 

Frequency multiplication is accomplished 
experimentally by individual frequency-doubler stages con- 
sisting of a full-wave linear rectifier followed by a 
bandpass filter tuned to the second harmonic of the applied 
signal. Because of the nonlinearity of this process, the 
output contains the product of signal and noise components. 
These additional terms reduce the output SNR. An analysis 
of the SNR reduction for the full-wave linear rectifier is 
difficult [251. On the other hand, the noise performance 
of a full— wave square— law device is well known [26] . Since 
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It has been shown experimentally that the results obtained 
using linear and square-law detectors are quite similar 
[6, p. ^32], it is reasonable to use the following 
analytical results obtained assuming a square-law frequency 
doubler . 



The SNR at the output of each doubler as a 
function of the input SNR is given by 



r. 

r = ^ 

o 2(l+2r^) 



For large this becomes 



(i^8) 




which is equivalent to a 6-dB loss in SNR per stage. This 
loss factor is only an estimate since it assumes square-law 
processing of an unmodulated carrier plus noise, whereas 
the doppler signal is frequency modulated. Further, an 
increase in SNR may result from bandpass filtering the 
output of the last frequency-doubler stage. In general, 
this improvement by filtering arises because the usable 
signal power tends to be concentrated in the filter pass- 
band while in most cases of practical interest the prominent 
noise terms after frequency multiplication have frequencies 
which are outside the filter passband. 

It is of interest now to consider the effect 
on the measurement of azimuth angle of the loss in SNR 
due to frequency multiplication. Based on the results 
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given by Skolnik [6, p. 462-478] for measurement accuracy 
under high SNR conditions, the angular error e(6) v/hlch 
arises is of the form 



(S/N) ^ 

Where 0g Is the half-power width of the output 
pulse from the MFD DP system with no frequency multiplica- 
tion and K is a constant. 

Now consider a system consisting of a single 
frequency doubler preceding the bandpass filter and assume 
that the 6-dB SNR loss caused by frequency multiplication 
is reduced to 3 dB by the filtering explained previously. 
Then at the system output the angular error e'(6) is 
given by 



s' (6) 



% 



( S ' /N ' ) ' 



Since S/N = 2(S'/N')j then 



(51) 



e'(0) 



(S/N)'^ 



(52) 



which means that if the new angular error e' is to be equal 
to or less than the angular error e which exists without 
frequency multiplication, then 



®B - ^ 

B /2 
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For a system containing several frequency - 
doubling stages, the 6-dB loss in SNR per stage requires 
that the beamwidth be halved with each stage to preserve 
the same accuracy. Such improvement is not possible as 
shown in Fig. 15. 

The preceding analysis estimates the effect 
of noise on a system using frequency multiplication. The 
disadvantages which arise from frequency multiplication may 
be balanced by other factors such as: 

(1) The accuracy obtainable by a human operator is 
influenced by the width of the pulse which must be "split" 
by a cursor or other bearing indication device particularly 
at large values of SNR. An operator can position a cursor 
more accurately on a narrower pulse thus counteracting 
some loss in theoretical accuracy due to multiplication. 

(2) Integration of a number of sweeps may increase the 
SNR without increasing the width of the response main lobe. 

(3) Additional pulse compression achieved at higher 
values of 6 which are obtained by frequency multiplication 
should offset some of the loss in SNR caused by the 
multiplication process. 

The next section gives a partly qualitative and 
partly quantitative confirmation of the analyses presented 
in this section. 

b. Experimental Results 

This section presents the results of a number 
of experiments in which signal plus additive bandllmited 
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noise is applied to an experimental MFD DF system (Fig. 

38 ) using frequency multiplication by a factor of four. 

Two photographs are generally shown for each input. The 
first photograph shows two cycles of the output of the 
final bandpass filter (DF system output). The second 
photograph shows the system output after envelope detection 
by a diode followed by an RC lowpass filter having a time 
constant of ^.7 ms. 

In the first experiment, bandllmlted gausslan 
noise having the spectral shape shown in Fig. 39-a is 
added to the FM signal. The value of 3 is 25. The signal 
plus noise is frequency multiplied by a factor of four to 
produce a signal with a value of 3 of 100. Mixing and then 
filtering with a 4-kHz bandpass filter complete the system 
operations. The results are summarized in Table IV and 
Fig. ^5 for a number of Inputs of differing values of SNR. 
The reduction in SNR by frequency multiplication is seen 
by comparing these results with the first non-multipllcatlve 
case (3 = 25) given by Table I and Fig. 42. Because of the 
interaction of signal with noise in the frequency- 
multiplication process, measurement of the output signal 
power (with no input noise), output noise power (with no 
input signal), and output signal-plus-noise power separ- 
ately provides no valid indication of the extent to which 
the signal Interacts with the noise. Therefore, systems 
which use frequency multiplication cannot be compared 
directly with those that do not on the basis of output SNR. 
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TABLE IV 



EFFECT OF NOISE ON THE PERFORMANCE OF 
AN EXPERIMENTAL MFD DF SYSTEM 



Basic Parameters : 

(a) B (before frequency multiplication) = 25 

(b) B (after multiplication In frequency by a factor 

of four) = 100 

(c) ^-kHz bandpass filter (TL iJD8-A) 



Figure 


INPUT 




OUTPUT 


Signal 

(dB) 


Noise 

(dB) 


SNR 


Signal 

(dB) 


Signal plus 
Noise (dB) 


^5-a,b 


8.0 


5.0 


3.0 


-20.0 


-16.7 


^5-c,d 


8.0 


-2.0 


10.0 


-20.0 


-17.2 


^5-e,f 


8.0 


-12.0 


20.0 


-20.0 


-19.0 
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(a) 




(b) 



FIG. 45. PHOTOGRAPHS OF THE RESPONSE OF THE EXPERIMENTAL 
MFD DF SYSTEM FOR VARIOUS VALUES OF INPUT SNR 
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(c) 




(d) 
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(e) 




(f) 
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However, a visual inspection of the photographs for the two 
cases gives a good indication of the deterioration produced 
by frequency multiplication. For example. Fig. i|2-e,f, the 
3-dB input SNR case, is visually comparable to the 10-dB 
input SNR case in Fig. ^5-c,d. This comparison indicates 
that the presentation obtained using a system without 
frequency multiplication is comparable to that obtained 
with a times-four frequency-multiplying system having a 
7-dB larger value of input SNR. Comparisons of other 
photographs generally confirm the 7-dB figure. 

The second experiment of interest is identical 
to the first except a 10-kHz bandpass filter is used 
Instead of the 4-kHz bandpass filter. The disadvantage 
of the increased bandwidth is evidenced by the results 
given in Table V and shown in Fig. which show on the 
average a 3-to 4-dB loss in the SNR of the system output 
relative to the results given in Table IV. 

In a third experiment, the tuned circuit 
which filters the wideband noise spectrum at the output 
of the noise generator, is replaced by a 10-kHz ceramic 
bandpass filter and the 4-kHz ceramic filter is used to 
select the desired frequency group. The results, summarized 
in Table VI and Fig. 47, are surprising because the noise 
tended to concentrate around the output pulses much more 
than when tuned LC circuits are used to limit the noise 
bandwidth. There is no obvious explanation for this 
phenomenon . 
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TABLE V 



EFFECT OF NOISE ON THE PERFORMANCE OF 
AN EXPERIMENTAL MFD DF SYSTEM 



Basic Parameters : 

(a) 3 (before frequency multiplication) = 25 

(b) 3 (after multiplication in frequency by a factor 

of four) = 100 

(c) 10-kHz bandpass filter (TL 10D9-20A) 



Figure 


INPUT 




OUTPUT 


Signal 

(dB) 


Noise 

(dB) 


SNR 

(dB) 


Signal 

(dB) 


Signal plus 
Noise (dB) 


46-a,b 


8.0 


5.0 


3.0 


-18.0 


-10.7 


46-c,d 


8.0 


-2.0 


10.0 


-18.0 


-12.7 


il6-e,f 


8.0 


-12.0 


20.0 


-18.0 


-15.9 
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(b) 



FIG, 46. PHOTOGRAPHS OF THE RESPONSE OF THE EXPERIMENTAL 
MFD DF SYSTEM FOR VARIOUS VALUES OF INPUT SNR 
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(c) 
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TABLE VI 



EFFECT OF NOISE ON THE PERFORMANCE OF 
AN EXPERIMENTAL MFD DF SYSTEM 

Basic Parameters : 

(a) 8 (before frequency multiplication) = 25 

(b) 8 (after frequency multiplication by a factor 

of four) = 100 

(c) 4-kHz bandpass filter (TL i|D8-A) 



(d) 


10-kHz bandpass noise 


filter 


(TL 10D9 


-20A) 






INPUT 




OUTPUT 


Figure 


Signal 

(dB) 


Noise 

(dB) 


SNR 

(dB) 


Signal 

(dB) 


Signal plus 
Noise (dB) 


4?-a,b 


8.0 


5.0 


3.0 


-20.0 


-18.0 


0 

1 


8.0 


-2 . 0 


10.0 


-20.0 


-19.0 


4?-e,f 


8.0 


-12.0 


20.0 


-20.0 


-19.7 



1^11 




(a) 




(b) 



FIG. 47. PHOTOGRAPHS OF THE RESPONSE OF THE EXPERIMENTAL 
MFD DF SYSTEM FOR VARIOUS VALUES OF INPUT SNR 
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(c) 




(d) 
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E. EXPERIMENT USING A CONVENTIONAL DOPPLER DF SYSTEM 



This section describes the equipment, procedures and 
results of an experiment In which a MFD DF system was 
assembled by interfacing the hardware constructed for the 
laboratory experiments with a conventional doppler DF set. 
The experimental system was tested both with controlled 
target transmitters and distant targets of opportunity. 
Despite the discrepancies between the theoretical and 
experimental results, the experiment shows the MFD DF 
system can be used with existing commiutated circular 
arrays to provide Information on signal AOA. 

1 . Description of Equipment 

An AN/TRD-15 DF set [271 is used to provide live 
signals to test the MFD DF principle. The particular set 
used consists of an array of 25 monopoles equally spaced 
around a circle of 150 ft diameter. A matched cable from 
each antenna is connected to a stator of an RF rotary 
coupler physically situated In a hut at the center of the 
array. The effective point of reception is made to travel 
around the array by rotating at a 50 Hz rate a pickup- 
winding commutator located on the rotor of the coupler. 

The commutator rotation simulates the rotation of a single 
antenna around the circle. 

The rotor output is applied to a superheterodyne 
receiver and converted to an IF of ^55 kHz. In the conven- 
Uonal set, the IF is processed and eventually sent to a 
display unit. For this experiment, a parallel connection 
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was made from the last limiter stage of the receiver to the 
MFD DF experimental hardware. This allows both DF systems 
to function simultaneously. 

The experimental hardware consists of frequency 
multipliers, a balanced mixer, a ceramic bandpass filter, 
and optionally, an envelope detector (see Fig. 48). 

By varying the local-oscillator frequency, the 
difference frequency in the mixer output is shifted at will 
and any portion of the FM spectrum can be placed in the 
filter passband. The filtered output can be viewed on an 
oscilloscope, either before or after envelope detection. 

2 . Results 

During the course of the experiment several signals 
in the HF region were examined. They included controlled 
target transmitters, as well as targets of opportunity 
such as WWV. Because circuitry to remove modulation from 
incoming signals was not used, primary emphasis was on 
unmodulated or slowly modulated carriers. 

The results presented in Table VII and Figs. 49-53 
confirm some of the important characteristics of the MFD DF 
technique including 

(1) Filtering of the appropriate side frequencies of 
the output of a commutated circular array which operates in 
an equivalent manner to a single rotating element yields an 
output having an envelope consisting of a periodic pulse. 
The relative time of occurrence of the pulse is linearly 
related to the AOA of the received signal. 
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PERFORMANCE OF AN EXPERIMENTAL MFD DF SYSTEM 
USED WITH A CONVENTIONAL DOPPLER DF SET 
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FIG. 48. BLOCK DIAGRAM OF THE EXPERIMENTAL MFD DF SYSTEM 
INTERFACED WITH A CONVENTIONAL DOPPLER DF SYSTEM 



(2) The 180 ° shift In apparent bearing which occurs 
when the group of in-phase frequency components below the 
carrier frequency is selected is demonstrated. 

(3) The relative time of occurrence of the pulse 
(bearing) is shown to be independent of transmitter fre- 
quency . 

The system output for signal 1 of Table VII is 
shown in Fig. ^9 for approximately three cycles of 
scanning. The second pulse is weaker than the first and 
third because of incidental amplitude modulation of the 
envelope at twice the rotor revolution rate. The oscillo- 
scope from which the photographs were taken was set to 
trigger at a north reference; hence by relating the distance 
between pulses to one revolution of the rotor, approximate 
bearings can be derived from the position of the pulse 
along the horizontal axis. For example, in Fig. HS-a. the 
pulse appears at a point slightly more than one-quarter 
of a cycle from the triggering point; thus the signal 
bearing is estimated to be about 100 degrees. Figure ^9-b 
is an expanded version (horizontal trace) of Fig. 49-a 
showing the structure of the pulse in more detail. 

Signal 2 of Table VII (Fig. 50) shows the system 
response (three cycles) to an unidentified target of oppor- 
tunity. The main pulses in Fig. 50-b are narrower than 
those in Fig. 50-a since multiplication by a factor of four 
is used in the former case and by a factor. of two in the 
latter . 
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(a) Response of experimental MFD DF system to signal #1. 
(Table VII). Three cycles of scan. 




(b) Response of experimental MFD DF system to signal (/I. 
(Table VII). Part of one scan. 



FIG. 49. PHOTOGRAPHS OF THE RESPONSE OF THE 
MFD-AN/TRD-15 DF SYSTEM 
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(a) Response of experimental MFD DF system to signal #2 
(Table VII). 




(b) Response of experimental MFD DF system to signal #2 
(Table VII), 



FIG. 50. PHOTOGRAPHS OF THE RESPONSE OF THE 
MFD-AN/TRD-15 DF SYSTEM 
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Figs. 51-53 are photographs of detected voltages. 
For example. Fig. 51 shows the envelope of the system out- 
put for an unidentified signal (number 3 of Table VII). 




FIG. 51. PHOTOGRAPH OF THE ENVELOPE OF THE RESPONSE 
OF THE MFD-AN/TRD-15 DF SYSTEM 

As discussed previously in Chapter II, the frequency 
components of the spectrum of a carrier frequency modulated 
by a sinusoid are symmetrical with respect to magnitude, 
though not phase, about the carrier frequency f ^ . Filter- 
ing the in-phase group of frequency components below f^ was 
shown to produce a pattern identical in shape to that obtained 
by selecting the corresponding group above f^. The only 
difference in the outputs is a time delay corresponding to 
a half period of rotation of the rotor, or in terms of AOA 
measurement, a bearing shift of l80 degrees. Despite the 
poor quality of the pulses. Fig. 52 illustrates this result. 
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(a) Detected response of experimental MFD DF system 
to signal #4 (Table VII). In-phase group above the 
carrier frequency. 




(b) Detected response of experimental MFD DF system 
to signal #4 (Table VII). In-phase group below the 
carrier frequency. 



FIG. 52. PHOTOGRAPHS OF THE ENVELOPE OF THE RESPONSE 
OF THE MFD-AN/TRD-15 DF SYSTEM 
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Fig. 52-a shows the transmitter bearing to be about 100 

degrees when the in-phase frequency components above f are 

c 

filtered, while Fig. 52-b shows the same signal at an 
apparent bearing of about 280 degrees when the in-phase 
group below f^ is filtered. 

Fig. 53 is intended to show that the presentation 
obtained with a target transmitter at a fixed bearing is 
unchanged as the frequency changes from 21.7^ MHz to 29.95 
MHz. 

Discrepancies between the results Just presented, 
the results obtained from laboratory experiments, and the 
computer results can in part be attributed to the following: 

(1) The inductive sampling of RF energy from 25 discrete 
antennas instead of a continuous output from a single 
rotating antenna. 

(2) Incidental FM caused by the limiter. 

(3) Non-ideal amplitude- and phase-vs-frequency charac- 
teristics of the RF and IF portions of the superheterodyne 
receiver. 

(^) Lack of knowledge of the position of the received 
signal in the passband of the receiver. The local 
oscillator was tuned until an optimal or near optimal MFD DF 
system output was observed on the oscilloscope. Use of a 
spectrum analyzer would have permitted more precise tuning. 
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(a) Detected response ot experimental MFD DF system 
to signal #4 (Table VII), Transmitter at 21.74 MHz. 




(b) Detected response of experimental MFD DF system 
to signal #5 (Table VII). Transmitter at 29.95 MHz. 



FIG. 53. PHOTOGRAPHS OF THE ENVELOPE OF THE RESPONSE 
OF THE MFD-AN/TRD-15 DF SYSTEM 
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Particular discrepancies which appear in Figs. ^9-53 are 



a. pulses narrower than anticipated. 

b. undesired secondary pulses 

c. irregular pulse shapes. 
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V. THE TWO SIGNAL CASE 



A common occurrence In the propagation of radio waves 
is the arrival of a signal at a receiver by two different 
paths. This phenomenon Is particularly troublesome in DF 
when the multipath signal arrives from two different 
azimuth angles. This chapter explores analytically the 
theoretical performance of a MFD DF system under multipath 
conditions. The effect of angle of elevation is briefly 
considered . 

A . THEORY 

When two signals of the same frequency with angles of 
arrival 0^ and 02 are received by an element (or equivalent 
array) oriented initially as shown In Pig. 5^ and rotated 
(or sampled) counterclockwise, the response of the system 
Is 



p(t) = C^cos[(jo^t + 3sln(oj^t - 0^)] 

+ C-ccc [ou t + B3ln(Lo t - 0.,) + 1,] (5'.' 

2 c r 2 ^1 

where is the phase difference between the Incoming 
signals when received at the center of rotation and and 
C-. are constants. 

cL 

To simplify (53) let 

Ci>C 2 , = 1 , b = 02 / 0 ^ , and 0^ = 0. 
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FIG, 54. PLAN VIEW OF A ROTATING ANTENNA ELEMENT 
SHOWING THE RECEPTION OF TWO LIKE -FREQUENCY 
SIMULTANEOUS SIGNAl.S 



Under these conditions ( 53 ) becomes 
p(t) = cos[(jo^t + Bsinoj^t] +b cos [w^t + + 6 sin(w^t - 62) ] (5^) 

Because of the symmetry in the spectrum of sinusoidal 
FM, no generality is lost by considering only the frequen- 
cies above the carrier. Following a development similar to 
the single-signal case in Chapter II, and shown in detail 
in Appendix B, two groups of in-phase harmonics, one from 
each term in (5^) 3 can be isolated and expressed in the 
form of a carrier and a modulating envelope as 
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p^(t) = F(t)cos[oj^t + <t(t)] 



(55) 



where F(t) Js the arrp3 itude-nodulaticn factor aiid C ‘ lo 
the phase- modulation term. 

Concerning MFD DF processing of mu'’ Ipati. signals it is 
of interest to determine whether the envelope F(t) at the 
output will have two main pulses corresponding in tim.e to 
the azimuth angles of the two incoming signals, assuming 
that the azimuth separation is sufficient for resolution. 
From (55) and Appendix B it appears that there are tvjo 
terms which may correspond to two separate pulses. The 
only uncertainty is the effect of the phase difference 

Figure 55 confirms the existence of two separate pulses 
in the case of two equal-am.plitude signals, one et 




FIG. 55. COMPUTED ENVELOPE OF THE RESPONSE OF AN IDEAL 
MFD DF SYSTEM TO TWO SIMULTANEOUS SIGNALS 
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l80 degrees of azimuth and the other at 270 degrees, ^^rith 
a relative carrier phase difference of zero. In this 
example, a value of S of 50 Is assumed as well as Ideal 
filtering of the In-phase group of frequencies (44th 
through 58 th harmonics of f^ above ■ 

The effect of the phase difference can be seen by vary- 
ing fi’om zero to 360 degrees. Calculated results show 
that causes the peak-to-dlp ratio to vary slowly from 
9.6 dB when = 253 degrees to 51.4 dB when Tpj = 70 degrees. 
However, variations in the location of the two maxima, 
corresponding to errors In bearing measurements are less 
than 1.6 degrees for all values of Figure 56 shows 

the "best" and "worst" cases for this example corresponding 
to the two extremes In peak-to-dlp ratio mentioned above. 

The MFD DF system response to two signals of the same 
frequency thus consists basically of the superposition of 
each individual response. The phase difference between 
the two arriving signals may affect the interaction of the 
two responses . 

The relative amplitudes of the two received signals 
can reasonably be expected to directly affect the relative 
amplitudes of the two pulses at the output of the MFD DF 
system. This Is exemplified in Fig. 57 which shows computed 
output envelopes from a system having a value of S of 50 
for the case of two signals separated 90 degrees in azimuth. 
Four different relative amplitudes are shown In Fig. 57 
ranging from b = 0.9 to b = 0.1. The error in the derived 
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(b) Phase difference is 253 degrees. 



,FIG. 56. CCMPUTED ENVELOPES OF THE RESPONSE OF AN IDEAL 
MFD DF SYSTEM TO TWO SIMULTANEOUS SIGNALS 
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o.2u 3 . g : 



(a) b - 0.9 




(b) b - 0.7 



JFIG, 57. COMPUTED ENVELOPES OF THE RESPONSE OF AN IDEAL 
MFD DF SYSTEM TO TWO SIMULTANEOUS SIGNALS 
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(c) b * 0.4 




(d) b * 0. 1 
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azimuth angle of the weaker of the two varies from 1.0 
degree (b=0.9) to 12.5 degrees (b=0.1). 



As the value of 3 Increases the system output pulses 
which Indicate the AOA of the received signals become nar- 
rower. This Is exemplified by Fig. 58 which shows the 
response produced by two signals 90 degrees separated In 
AOA when the value of 3 Is 200. 



FIG, 58. CCMPUTED ENVELOPE OF THE RESPONSE OF AN IDEAL 
MFD DF SYSTEM TO TWO SIMULTANEOUS SIGNALS 

B. RESOLUTION 

The theoretical resolution of a DF system Is a measure 
of the smallest difference In AOA between two arriving 
signals at which each can be separately distinguished. Fo 
the output of the MFD DF system, this corresponds to some 
measure of the dlstlngulshablllty between the two pulses 
corresponding to the AOA of the two received signals. The 




A 



A 



measure taken here is the smallest interval (expressed in 
degrees of azimuth) for which the dip between the two 
pulses is at least 3 dB below the smaller of the two. 

The most Important parameters affecting resolution are 
relative signal amplitude b, carrier phase difference , 
and modulation index 3- 

The effect of the relative amplitude b is to decrease 
the resolution as b decreases. Thus, as shown in Fig. 59 
for 3 = 50.0 and = 0 the resolution angle varies from 
60.7 degrees when the two signals are of equal amplitude, 
to 102 degrees when the amplitude of the weaker signal is 
0.1 that of the stronger (20-dB level difference). 

The effect of the carrier phase difference is seen 
in Fig. 60 obtained at a constant value of 3 = 50 for 
equal-amplitude signals, which shows the resolution angle 
varying from 55-5 degrees to 6 I .5 degrees as goes from 
0 to 27 t radians . 

Finally, the effect of 3 on resolution is presentea in 
Fig. 61 as the other two pertinent parameters b and are 
held constant. As might be expected from previous dis- 
cussions the resolution improves directly with 3. 

C. INFLUENCE OF ELEVATION ANGLE 

The effect of the angle of elevation a of the received 
signal on the FM signal produced by a rotating antenna 13 
can be taken into account by rewriting (5^) as 

p(t) = cos[oj t + 3cos a, sin(cj t ) ] +bcos [oj t+i|^, +3cos apSin(w t-6p) ] 

C J_X* C-L d. L ^ 

(56) 
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FIG, 59. GRAPH OF THE THEORETICAL R SOLUTION AND BEARING ERROR OF THE MFD DF SYSTEM 
VS THE RATIO OF THE AMPLITUDES OF TWO LIKE -FREQUENCY SIMULTANEOUS SIGNALS 




saaa^^p ‘uoT;n]Osaa .iexn§uy 



c. 60. GRAPH OF TH*^ THFORFIfCAI RFSOII'TTON TII‘ NFD IR^' YSI'Ol VS THF 
I'lFFERENGF , IN PHASE OF HIF IWo I.IKE-FRFQL’f NGY STmL’LI A.N OUS SKJNAES 



o 

o 

00 




I 

o 

00 



— J I L i- A I 

o o o o 

vo CM cn 

S33aSap ‘uofqn-josaa aB'[nSuv 



168 



:G. 61. GRAPH OF THE THEORETICAL RESOLUTION OF THE 
MFD DF SYSTEM VS THE MODULATION INDEX B 



The effect of elevation angle can be equated to a reduction 
in the modulation index 3. As long as and are 
approximately equal, the effect is a decrease in the effec- 
tive aperture of the antenna since a lower value of $ can 
be associated with a smaller radius of rotation. However 
If and differ significantly, there will be a consider- 
able difference in the overall bandwldths of the two FM 
signals at the output of the antenna. This means that a 
bandpass filter which selects the in-phase group of fre- 
quency components from the signal having the larger value 
of 3 (smaller a) will have fewer or none of the desired 
spectral components of the signal having the lower value 
of 3 (larger a) in its passband. The MFD DF system output 
under these circumstances de-emphasizes the effect of the 
signal having the largest elevation angle or would provide 
no indication of its presence. 

On the other hand, to attempt to bandpass filter the 
in-phase group of frequency components associated with the 
signal having the smaller value of 3 would necessarily 
admit undesired frequency components of the signal having 
the larger value of 3, so this procedure may not be 
acceptable . 
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VI. SUMMARY AND RECOMMENDATIONS 



A. SUMMARY 

The problem of determining the direction of arrival of 
a radio wave has been approached from a different point of 
view, namely, that of applying matched-filter processing 
to the output of a rotating antenna. 

Chapter II presents the principal results of this 
research which are the development of the theory and an 
indication of the performance of a matched-filter doppler 
DF system which is simple to realize and completely linear. 
When a radio wave impinges upon an antenna element omni- 
directional in the plane of propagation and which rotates 
in a circle (also in the plane of propagation) the output 
of the antenna is modulated sinusoidally in frequency due 
to the doppler effect. If this modulated signal is applied 
to a bandpass filter which allows only a particular group 
of frequency components to pass, the output is a periodic 
train of short duration pulses having no sidelobes. The 
time of occurrence of each output pulse relative to one 
rotation period is directly related to the AOA of the radio 
wave . 

This technique evolves as a practical compromise after 
recognizing the difficulty in realizing a filter matched 
to the entire spectrum of a signal frequency modulated by 
a sinusoid. The bandpass filter Isolates a group of 
frequency components which by themselves satisfy without 
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processing of any kind the condition of in-phase frequency 
components normally obtained with a matched filter. 

To Improve system accuracy it is desirable that the 
duration of the pulse used to measure AOA be small. Pulse 
duration is inversely related to the value of the modulation 
index 6 of the FM signal. In turn, at any given frequency, 
the value of g is directly proportional to the radius of 
the circle of rotation (or of an equivalent circular array). 
It appears, therefore, that a reduced pulse duration results 
only from an Increase in the size of the receiving aperture. 
However, it is shown that the technique of frequency multi- 
plication can be used to Increase the value of g without 
increasing the aperture size. 

To apply signal-processing techniques, such as matched 
filtering, to obtain spatial information, such as direction 
of arrival, the receiving antenna must convert spatial 
Information to temporal form. As shown in Chapter III this 
is exactly what occurs when a linear antenna (or array) is 
rotated about a central point. This chapter also Includes 
a Fourier transform representation of an antenna and 
reveivlng system. 

The physical realizability of the MFD DF technique. 
Including the frequency-multiplication scheme, is demon- 
strated by presenting results obtained with a hardware 
simulation system. A critical component is found to be the 
bandpass filter since Imperfect filtering of the desired 
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frequency components is found to produce undesireable 
results such as generation of double maxima in the output. 

The effect of bandlimlted gaussian noise on the MFD DF 
system is considered. When processing without frequency 
multiplication the system Improves SNR due to the 
difference in spectral distribution of the noise power and 
the signal power. There is also an increase in the ratio 
of peak signal to average-noise power because the output 
signal power is concentrated in a fraction of the cycle 
(pulse compression) while the noise power is distributed 
over the entire cycle. When the processing of the signal 
includes frequency multiplication (a nonlinear step), the 
SNR is reduced by about 6 dB each time the frequency is 
doubled. This loss in SNR tends to degrade system accuracy 
whenever the input SNR is less than about 20 dB. There- 
fore, frequency multiplication may be a useful option in 
system operation only if the input SNR exceeds about 20 dB. 

An experiment in which the experimental MFD DF cir- 
cuitry is Interfaced with a conventional doppler DF set 
confirms under field-operating conditions the basic 
characteristics of the MFD DF technique. 

The theoretical response of the system to two signals 
of the same frequency at different angles of arrival con- 
sists basically of the superposition of the two individual 
responses . 
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B. RECOMMENDATIONS 



1. To further evaluate the merits and limitations of 
the MFD DF technique an operational system should be con- 
structed and tested. 

2. A tactical DF system based on the MFD technique 
may be possible by Installing an antenna in the end of a 
helicopter rotor. 

3 . A study should be made of the spectrum of various 
types of frequency or phase modulation with the object 

of finding signal spectra which are amenable to matched- 
filter processing using simple conventional filters. 

4. Other methods of direction finding should be 
studied to determine if various means of increasing signal 
bandwidth can be used advantageously in conjunction with 
these DF techniques. 
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APPENDIX A 



COMPUTATION OP THE RESPONSE OF AN IDEAL 
MATCHED-FILTER DOPPLER DF SYSTEM 



The response of an ideal doppler DF array of radius i 
to an unmodulated signal at frequency f^ and azimuth 6^ 
is given by (1) as 



s(t) = Ccos[u) t + sln(o) t - 6 )] 

'• c A r o 



(A-1) 



Expanding the above in a Fourier series yields (when 
C = 1), 



s(t) = [J (3)costo t+ E J, (3)cos[(o) +kto )t-k6 ] 
o c ^ c r o 

+ E (-1)^ Ji^( e)cos [( w - kw )t+ke ] (6) 

k=l ^ c r o 



Because Jj^(3)-^o as k>>3 the above series may be trun- 
cated at some integer n^ for which (3) has decreased to 
some small factor of the largest Jj^(3). The effect of an 
idealized bandpass filter is simulated by omitting all 
terms in the series except those which correspond to the 
frequencies which are passed by the filter. This yields 



s^(t) = E Jj^( 3)cos [ ( oj + ka)^)t - k0^] (A-2) 

k=n 

o 

To study the envelope E(t) of (A-l), a conversion is 
made, using trigonometric identities, as follows 
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(A-3) 



s,(t) = AT(t)coso) t - B,(t)slnto t 
i 1 cl c 



where 



A, (t) = E J (3) [cos k(,o) t - 0 )] 

k=n^ 

B, (t) = E J, ( 3) [sin k(oj t - 0 )] , 

^ k=n r o 

o 



Finally , 



s^(t) = E(t)cos[co^t + <J)(t)] 



where 



E(t) = 



Z[A^(t)] 



^ + [B^(t)]^ 



and 



cf)(t) = tan 



_1 B^(t) 



A^(t) 



(A-i^) 

(A-5) 

( A-6 ) 

(A-7) 

(A-8) 



Values of E(t) are obtained with the aid of a digital 
computer. 
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APPENDIX B 



COMPUTATION OF IDEAL MFD DF SYSTEM 
TWO-SIGNAL RESPONSE 



The response of an Ideal doppler array to two signals 
at the same frequency but different azimuth angles is given 
by 

p(t) = C^cos [co^t + 3 sln(o 3 ^t - 9j^)] 

+ C 2 C 0 S [o)^t + Bsln(o3^t - 62 ) + 4^^] (53) 

This expression can be simplified by assuming 
Cf = 1 , C 2 <C ^5 b = C^/C^, and 9^ = 0 . 

Under these conditions, 

p ( t ) = cos + 6 slnw^t ] + bcos[o 3 ^t + 3 sln( co^t - 92 ) ] (B- la) 

from which 

p(t)= cos[(jj^t + 3 slnw^t ] +bcos4^^cos[u)^t + 3 sin( co^t - 92 ) ] 

-bsln4^^sln[a)^t + 3 slnCw^t - 92 ) ] . (B-lb) 
Expanding (B-lb) in a Fourier series yields 
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n., 



n. 



1 1 

p(t) = J (B)cosco t+ E J, (6)cos(o) +kto )t+ Z J (3)cos(o) -kw )t 
o ^ k=l ^ ^ ^ k=l ^ c r 



’^1 

+ bcos'/^-| J (3)cosco t+ Z J,(B)cos[(co +k(jo )t-k6„] 
_L \ O C ^ 1 ^ 0 27 ^ 



n- 



k 



+ Z (-1) J, (3)cos[(o) -kta )t+ k0„] 
k=l ^ c r ^ 



r 

-bsln'l'. J (3)slnto t+ Z J, (3)sin[(co +kco )t- k6„] 

_L', O Ci_-iK C]T* d. 

I k=l 






,k 



+ Z (-1) J^(3)sin[ (oj^-kco^)t + ke^] 

k=l 



J 

(B-2) 



The effect of a bandpass filter Is simulated as in 
Appendix A by omitting all terms in the series except those 
which correspond to the frequencies which are passed by the 
filter. This yields 

*^1 

p (t) = Z J, (3)cos(oj + kco )t 
1 k=n ^ c r 

o 

'^l 

+ bcosijj, Z Jj^ ( 3 ) cos [ ( w^+kco^) t - kS^ ] 

k=n ^ 

o 

- bslnijj. Z J, ( 3)sln [ (to^+kco^)t - ke^] (B-3) 

k=n 

o 

To study the envelope F(t) of p^(t) (B-3) is put into a 
more useful form as follows 



p^(t) = A^it) cosixi^t - B 2 (t)sinto^t 



(B-^) 
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where 



A 2 (t) 


= E ( 6 ) { cos (kto t ) + bcos [k( ca^t - 02 ) + (B-5 ) 

k=n ^ ^ 

o 


B 2 (t) 


’^1 

= E ( 6 ) {sin (kto^t ) + bsin [k ( ui^t - 62 ) + 4 ^ 2 ^^ (B- 6 ) 

k=n 

0 



Finally 



P3_(t) 


= F(t)cos[co^t + (B-7) 



where 



F(t) = 


/[A2(t)]2 + [B2(t)]2 



and 



^(t) = 


B2(t) 

arctan ^ (B 9) 
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APPENDIX C 



MODIFICATION OF MFD DF TO PROCESS MODULATED SIGNALS 

The following is a summary of known means of removing 
modulation from a received carrier. 

A. AMPLITUDE MODULATION 

A limiter is sufficient to remove envelope variations 
of an AM signal 

B. ANGLE MODULATION 

For FM or PM signals the process consists of using an 
auxiliary receiver at the center of the array. The IF out- 
put of this receiver is 



where A(t) is the carrier angle modulation generated by the 
transmitter. 

The receiver associated with the doppler DF system 
yields an IF output given by 



s (t) = K cos[oj.,t + A(t)] 

3 . 3. _L 



(c-i) 



s,(t) = K, cos[w„t+A(t)+3sln(ojt-6 )] . 

D D ^ r o 



(C-2) 



If s (t) and s, (t) are applied to a mixer, the 

Si D 



difference-frequency output signal is given by 



s^(t) = cos [ )t + 3sln(oj^t - e^) ] 



(C-3) 
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which contains only the doppler modulation. The signal 
s^(t) is then processed for DF purposes in the usual 
manner. This technique of removing carrier angle modula- 
tion is used in conventional doppler DF systems [27] . 
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